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PREFACE 


This manual, describing the computerized coaxial injection combustion 
model, was prepared in support of the continuing JANNAF effort to 
develop systematic performance prediction techniques. 

In 1265, the Interagency Chemical Rocket Propulsion Group (ICRPG) Work- 
ing Group was formed for the purpose of improving and recommending method- 
ology suited to eventual adoption as national standards for the analytical 
and experimental evaluation of the performance of liquid propellant rocket 
engines. By 1968, the Working Group (made up of members from government, 
industry, and academia) had: 

Developed a physical model of rocket engine thrust chamber performance 

Selected computer programs to perform the mathematical calculations 
required by the physical model 

Developed recommended practices for test measurements 
Developed a model for uncertainty in measurements 

Documented the effort in three procedures manuals (CPIA No. 178, 179, 
and 180) and several computer program manuals. 

In 1968, the ICRPG was reincorporated as part of the Joint Army-Navy-NASA- 
Air Force (JANNAF) Interagency Propulsion Committee. The major JANNAF 
achievement to that time was the publication of standard Thermochemical 
Tables for rocket exhaust products. The ICRPG Performance Standardization 
Working Group became the JANNAF Rocket Engine Performance Working Group. 
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Other JANNAF Working Groups cover such areas as Combustion Instability 
(originally part of the ICRPG), and Air-Breathing Propulsion and Environ- 
mental Effects. Each Working Group has a four-person steering committee 
(each Government agency being represented) , a program manager to coordinate 
the Group's efforts, members from Government agencies, and participants 
from outside the Government. 


Since the reinstitution of the Rocket Engine Performance Working Group in 
1968, this Working Group has: 

Extended the methodology from the thrust chamber to the entire engine 

Developed a detailed injector analysis procedure to replace the earlier 
ICRPG method 

Developed a rigorous step-by-step analysis procedure and a simplified 
procedure using efficiency definitions 

Replaced the approximate boundary layer model with a more rigorous 
model 

Established new overall procedures and documentation consisting of a 
Performance Prediction and Evaluation Manual, a User's Guide based 
upon experience pertaining to the manual and recently a CPI A publi- 
cation (245) dealing with JANNAF Rocket Engine Performance Test Data 
Acquisition and Interpretation. 

Continued to update and improve all methods and procedures. 


The documentation of the coaxial injection combustion model contained in 
this manual is indicative of constant updating and improvement to JANNAF 
performance prediction procedures. Specifically, this report describes the 


use of a reference computer program developed for the rigorous analysis 
of rocket thrust chambers with coaxial propellant injection. An earlier 
version of the model described herein was referenced in CPIA Publication 
245 (page 13. 2B) as the CSS model, a "coaxial element model that replaces 
LISP and 3DC for coaxial elements". This report describes an improved 
computer program which supersedes CSS. The improved model is named CICM. 

This report has been prepared in fulfillment of contract NAS8-29664 from 
the National Aeronautics and Space Administration. The effort was com- 
pleted during the period from 2 May 1973 to 15 April 1974. Mr. K W. Gross 
of the NASA Marshall Space Flight Center was the Technical Monitor. The 
Rocketdyne Program Manager was Mr. L. P. Combs, initially, and later 

Mi. J. Friedman. Dr. Robert D. Sutton served as the Rocketdyne Project 
Engineer. 
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INTRODUCTION 


The performance of a thrust chamber depends greatly on the combustion 
and injection processes including atomization, evaporation, and mixing. 
The individual processes are highly complex, especially for coaxial 
liquid -gas streams. Over a period of several years, Rocketdyne has 
developed an analytical model to simulate these processes; the model 
has been used extensively in the Space Shuttle engine development effort. 

The JANNAF Performance Standardization Working Group has directed the 
development of reference computer programs for evaluation of liquid 
rocket engine performance. Current capability includes the ability to 
simulate the behavior of various injection elements such as doublets, 
triplets, pentads, showerheads, etc., for liquid-liquid propellants 
with the JANNAF Distributed Energy Release (DER) computer program. 
However, the ability to simulate coaxial liquid-gas jet injection and 
combustion was needed in the JANNAF methodology. 

The objective of this technical effort was to provide NASA/MSFC with an 
operational FORTRAN IV computer program and associated documentation 
applicable to analysis of coaxial injection and combustion of liquid-gas 
propellants in the Space Shuttle engine. In addition, the computer 
program was to be structured to fit into the standardized JANNAF evalu- 
ation procedure, so that other engine systems using coaxial injection 
could be analyzed. These objective., have been met. 
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The effort was divided into two tasks. The first concerned improvements 
vo the existing Rocketdyne model: program modularization, improvement 
of the numerical analysis, modification of program tables, changes to 
the input and output format, and inclusion of punched card output compat- 
ible with the JANNAF DER program input requirements. The second task 
involved documentation of the model. 

This report has been prepared to provide sufficient information to allow 
MSFC to adequately use the model. This report includes overall descrip- 
tions of the equations, the overall program and its subroutines (includ- 
ing flow charts that emphasize the interactions of the subroutines rather 
than the detail of their internal structure), the program input and out- 
put, internal checks, guidelines, and error analysis. 

SCOPE AND LIMITATIONS OF THE 
COMPUTER PROGRAM (ClCM) 

This report describes a very comprehensive and complicated computer pro- 
gram to predict the combustion within rocket thrust chambers of gas/liquid 
propellants injected from coaxial elements. The model is capable, when 
used with applicable intra-element mixing data, of predicting the perform- 
ance of any size of concentric coaxial element using any propellant 
combination. 

The program is designed to use the injector and chamber configuration, the 
propellant and the operating conditions as the input. In a single run, the 
program will calculate the state of flow conditions within each element's 
"cup" (volume formed by recessing the oxidizer propellant delivery post 
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into the fuel sleeve). The calculated flow conditions include the 
important fuel side "cup' ; pressure drop. This pressure loss is 
defined to be the difference between the pressure in the fuel 
annulus jjap compared to that at the injector face. The program 
then automatically stores and uses all data as input to the chamber 
calculation sequence. If the injector lias more than one element de- 
sign (or zones representing a group of elements each having similar 
inlet conditions) the model performs automated repetitive analysis 
until the element (zone) with the longest predicted jet length is 
located. At .his point spray gas information from all of the elements, 
or groups of elements, are internally input into an auxiliary program 
that organizes the data, in terms of punched card output, for sub- 
sequent analyses of the completion of combustion, etc. via other 
computer programs ( J AXNAF DER.et scq.). The zones of element inlet 
conditions (i.e., injector feed maldistribution) must of course be part 
of the input to DER. 


This rocket engine combustion model is unique because it calculates both 
the rate of atomization of the injected liquid jet, resulting from the 
sheer force between the jet and the surrounding gas, and the axially 
varying mean droplet size produced by the atomization. Thus, it does 
not require experimentally determined correlations for the droplet size 
distribution, which are required in other models. 


Integration of the computer program into the DER methodology was formulated 
after considering many alternative methods for handling intra-element mixing. 
It was decided that the most accurate way to compute the total effect of 
this phenomenon is to divide the spray and gas flows for each coaxial element 
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zone into multiple mixture ratio sub-zones. The manner in which the 
element zones are further subdivided to simulate the intra-element 
mixing loss is determined from cold flow measurements. These cold 
flow measurements relate the element geometry and flow condition to its 
mixing efficiency. Such information may be input to the CICM program 
in terms of mass fraction as a function of the total fuel and oxidizer 
flowrate for each element zone. This analysis is used in an auxiliary 
program that interfaces CICM with the streamtube portion of DER. 

The streamtube portion (STC) of DER must, of course, be provided with 
more information than the punched card output from CICM provides. In 
essence, STC contains multiple concentric streamtubes representing each 
zone and the further breakdown of each zone into additional concentric 
streamtubes to account for intra-element mixing efficiency. Although 
mass fraction as a function of the fuel and oxidizer flowrates for each 
element zone is calculated, the user must decide, when inputing DER, which 
spatial concentric streamtube to use to represent each zone, and further, 
what mixture ratio profile to assign to the additional concentric stream- 
tubes within each overall zonal streamtubes. These same decisions are 
required for analysis of other element types when using the JANNAF DER 
(STC) program and, therefore, they do not represent additional complexity. 

The computer program also has capability for bypassing DER (STC) entirely 
and continuing the spray/gas combustion computations for single streamtubes 
to the nozzle throat. In such a case the area of the streamtube varies as a 
constant proportion of the total cross-sectional area, which is usually based 
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on tne iatio of element flowrate to total flowrate. Unlike DER (STC) , in 
constant area sections of the chamber the strenmtube area does not change, 
rather it retains a constant area. This simplification occurs because only 
one streamtube is being considered, whereas in STC many streamtubes are 

considered and adjustments are made to their individual areas so that they 
sum to the chamber area. 


All physical properties in the program are supplied by generalized property 
table subroutines for all droplet liquid and vapor, combustion gas, non- 
burning gases and droplet film properties. The program utilizes an advanced 
droplet vaporization and heating model which includes real gas effects re- 
garding vapor-liquid equilibrium and solubility of external gases into the 
droplet. To describe these non-ideal gas effects, the Redlich-Kwong equa- 
tion of state and the mixing rules of Chueh and Prausnitz* have been utilized. 
Rocketdyne has developed separate programs to calculate the non-ideal effects 
required in the CICM program. Although the program has been generalized to 
accept any propellant combination, non-ideal properties have been supplied 
only for the LiyGHj propellant combination. Additional effort would be 
required to supply properties for other propellants. 

During the current effort, the non-ideal physical properties for the LcyGH, 
propellant system were evaluated to ensure that the program is adequate*/ 
computations to at least 5000 psi. Non-ideal effects (where applicable) of 
temperatures from 100°R to 10.000°R have also been included. 


Chueh, P. L. , and J, M, Prausnitz. 
Liquid Equilibria," Industrial and 
1968, pp. 34-52. 


^Calculation of High-Pressure Vapor- 
Engineering Chemistry . Vol. 60, 
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As noted earlier, the program contains an advanced droplet vaporization 
and heating model (similar to, but an improvement on, that contained in 
the current DER/STC droplet heating program). This model was developed 
to permit analysis of both subcritical and supercritical temperature and 
pressure conditions. The model predicts continuous variation of burning 
rate with pressure. It computes a "wet bulb" temperature for subcritical 
pressures while allowing the droplet to continue heating through and past 
the critical temperature for supercritical pressures. 

The portion of the program that deals with the "cup" region (that volume 
created by recessing the oxidizer post within the fuel sleeve upstream 
of the injector face) permits analysis of both non-burning or ignited gas 
flows. Ignition of the injected gas (usually fuel) and atomized and 
vaporized oxidizer (the liquid jet) is believed to occur primarily as a 
result of recirculation of hot gas which is promoted by a highly flared 
fuel sleeve (such as on the J-2 and J-2S). With such a flare, an igni- 
tion front is established across the gas flow path at an angle determined 
by the cup gas propellant flow and flame speeds. Ignition is considered 
to occur at the beginning of the flare on the fuel sleeve and propagates 
downstream toward the liquid jet. For a non-flared cup, a similar analysis 
is used for ignition of non-burning gases as they enter into the chamber. 
Ignition occurs by local recirculation of hot gas around or between elements. 
Determination of the chamber gas flame speed is not exact; generally a value 
of 600 feet/second is recommended. 
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Other ignition mechanisms are possible but no high performance coaxial 
injected engine is known to have been (or is ever likely to be) built 
which has a gas injection velocity low enough to allow a turbulent flame 
to propagate upstream into the cup and maintain combustion there. 
Similarly, no engine is known to have been built or designed with gas 
injection temperatures high enough to cause auto-ignition to occur. 

Even in the case of the SSME, tests have been made with hydrogen-rich 
gas injected at 2000°R (considerably hotter than that planned for the 
actual engine) without indications of cup ignition. Ignition is un- 
likely because the induction time for auto-ignition at these extreme 
conditions is some 30 times greater than cup gas residence time. 


Additionally, the computer program has been modified from previous 
versions to improve the calculation of the "cup" pressure drop. With 
this modification, the pressure at the downstream end of the fuel gap 
annulus is predicted rather than at the propellant contact, point down- 
stream of the oxidizer post. This change is of significant importance 
to the accurate prediction of the overall fuel pressure drop, fuel mani- 
fold to injector face. The program will accurately predict the pressure 
recovery or loss (if any) in the flow from the annulus to the point of 
propellant contact, as well as the pressure drop from that point to the 
injector face. Three different methods are provided to compute the pres- 
sure differences between the point of contact and the fuel annulus. 

These three methods are discussed in detail under the General Program 
Outline section of this report. Of these three methods, one is recommended, 
but the others may be used with easily made changes to certain atomization 
rate and dropsize diameter correlation inputs which depend on the method 
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selected. Flame speed effects, whether in the cup or in the chamber, 
are relevant only with respect to use of the third (presently recom- 
mended) method. 

The primary limitation of the program is its essentially one-dimensional 
nature. Within any streamtube at a given axial location, the droplet 
spray and gas flow are considered uniform within the cross-section. 

Also, the program does not directly consider secondary droplet breakup. 
(Only the initial atomization rate and mean dropsize variation as a 
result of vaporization between axial steps are calculated.) However, 
neither of these restrictions seem to have any significant influence on 
the ability to model observed behavior. 

Non-uniformities are handled by the intra-element mixing technique des- 
cribed earlier. For the cup region, the method of predicting the atomi- 
zation rate and initial dropsize is believed to account implicitly for 
what secondary droplet breakup occurs because the rates and dropsizes 
are adjusted to obtain the best values of model parameters based on 
experimental information. The same statement can be made concerning 
the chamber, although for this region computations indicate secondary 
breakup is far less likely. In the chamber, the droplets are larger ini- 
tially because the flow is not constrained by the cup fuel sleeve walls. 
However, the computations indicate that the droplets rapidly vaporize 
in the surrounding hot combusting gas and are rapidly accelerated by 
droplet drag to a critical relative speed where further breakup (beyond 
the initial atomization) does not occur. 
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An additional restriction of the program is its inability to analyze 
coaxial injectors incorporating propellant swirl. 
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GENERAL PROGRAM OUTLINE 


Over a period oi several years, an integrated method for analyzing bipropel- 
lant liquid spray combustion has been developed and applied to steady-state 
wall heating and performance analyses (Ref. 1, 2, and 3). The approach was 
developed at Rocketdyne under a scries of contracts supported by the Air Force 
and by NASA and was guided by JANNAF Performance Working Group recommendations 
and requirements. For injection elements other than coaxial elements, this 
method is based on initializing the combustion field for the entire combustor 
(or a representative geometric segment of it) at a plane located a short 
distance downstream of the injector by summing the spray flux contributions 
from individual injection elements to each of a large number of flow-net mesh 
points. Individual element behavior is described analytically by empirical 
correlations of data from single-element, cold-flow experiments. Subsequent 
to that, combustion is described in a rapid combustion zone (if strong trans- 
verse gradients are produced) followed by a streamtube combustion zone, as 
shown in Fig. 1. The computer program (Ref. 3) necessary to analyze this 
fiOw field is schematically illustrated in Fig. 2. 

To date, this approach has been primarily applied to liquid-liquid systems. 

While cold-flow characterization of gas-liquid injection provides valid empirical 
correlations for many injector types, this approach may not be sufficiently 
accurate for injector types that exhibit strong coupling between the atomiza- 
tion processes and the spray combustion or vaporization processes. This is 
particularly true in regard to the dropsize distribution and location at which 
spray droplets are formed. Coaxial jet injection is subject to such coupling 
in at least two ways: (1) gas-liquid interaction in any elemental recesses 
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Figure 2 


Schematic of DER Computer Program for Performance 
Analysis 




te.g., such as recessing the liquid propellant injection post) accelerates 
the initial atomization, produces finer sprays, and increases injection 
pressure losses and (2) completion of jet atomization is accelerated and 
finer sprays are produced by the buildup of axial combustion gas velocity. 
Therefore, it was anticipated that some form of coupled atomization- 
combustion analysis of individual elements would be required for initializing 
the more global chamber combustion analyses. This is schematically implied 
in Fig. 3, The jet atomization may extend considerably into the combustion 
chamber, interacting strongly with (and producing) the surrounding 
combustion flow field. Thus, a jet atomization-combustion zone would 
replace the injection-atomization and rapid combustion zones of Fig. 1 . 

This requires an addition or alteration to the present computer program 
of Ref. 3, such as shown in Fig. 4, 


The Rocketdyne-developed coaxial injection combustion model (CICM) predicts 
the atomization, mixing, and (if present) combustion within the coaxial 

essed cup as well as the jet atomization and combustion within the 
combustion chamber. The basic program analyzes a single coaxial element 
or .element zone.' composed of similar elements which can be considered 
to be a single element. If manifold feed maldistribution is present (or 
if some of the elements differ in design) an internal multiple analysis 
is performed for each element zone that represents a different operating 
condition. Each element zone is assigned its proper proportion of the chamber 
are. and propellant flow. These element zones may be further divided 
into subzones to include the effects of intra-element mixing losses. 
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Proposed Restructuring of DER Program to Include Gas/ 
Liquid Coaxial Injector 


13 




ANALYTICAL FEATURES OF SINGLE-ELEMENT, COAXIAL 
INJECTION COMBUSTION 


As in most combustion analyses, inputs are required regarding the propel- 
lants' physical and thermochemical properties, equilibrium combustion pro- 
ducts, and injector, chamber and nozzle designs. The analysis of the 
injector element, as far as prediction of the vaporization efficiency is 
concerned, may proceed from a flow field formulation such as that depicted 
in Fig. 5. The element shown here is flared such as on a J-2 engine. (Most 
elements are not flared.) Only three axisymmetric flow fields need be con- 
sidered, as shown: (1) the liquid jet, (2) the spray/gas -burning (or non- 

burning) flow field surrounding the jet, and (3) bleed flow through the 
injector face, "Rigimesh flow", surrounding the gas/spray flow field and 
mixing with it. The flow within each of these fields is taken to be quasi- 
one-dimens ional (i.e., the radial mass flux concentration and pressure 
gradients are assumed to be insignificant). 

The actual analysis begins with initial contact between gas and liquid jet 
(Fig. 5). This contact may, for some injectors, occur in the cup formed by 
recessing the liquid oxidizer injection post. The conditions of this initial 
contact depend on the method chosen to describe the fuel flow from the fuel 
annulus gap to the liquid jet. In this region, the constrained high- 
velocity gas stream begins the initial stripping and atomization of the 
jet. Small drops are formed and, depending on local flow and geometric 
conditions (i.e., gas temperature, velocity, oxidizer vapor concentra- 
tion, pressure and flame speeds), the propellants may initially ignite 
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Figure 5. Conceptual Model of Uniformly Flowing Coaxial 
Injector Element 


withm the cup. Whether or not this occurs depends either on hot gas 
recirculation in the flare (if the element is flared) or on kinetic 
ignition delay times. Even if burning does not occur in the cup region, the 
recess of coaxial injectors has a significant effect on subsequent 
atomization and combustion in the chamber, thereby introlling the per- 
formance efficiency. In addition, significant cup burning often radically 
changes injection delta pressures and injection temperatures and must be 
accounted for in initial design. 

The program is based on the use of conservation equations for both the 
liquid jet and combustion gas/spray flow in the cup region and chamber 
flows. They include spray droplet atomization, heating, burning, and 
droplet drag, iieat transfer to the walls, injector facj, and liquid jet 
are neglected. In the chamber, uniform radial pressure is assumed at each 
axial location and the sum of the areas of the liquid jet, combustion/ 
spray field and "Rigimesh flow" must fill the portion of the chamber area 
allotted to the element. In the absence of "Rigimesh flow", the 
combustor flow field emanating from an element is allowed to expand at 
constant pressure until the flow fills the "chamber". This is not an 
assumption bur represents the limit of axial pressure variation as the 
ex .ernal (Rigimesh) flow is reduced to zero. These equations are described 
more fully later in this report. 
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Equations are also included which describe liquid jet stripping rates and 
drop-size distribution. This is the only known model which calculates 
stripping rates, atomization, and combustion of liquid jets in gas-liquid 
coaxial injectors. The current model represents a significant advance 
over the first attempt at modeling coaxial injectors (Ref. 4). The model 
does not require input data regarding dropsize information because the 
distribution is calculated as a function of flow field conditions and jet 
axial position. The controlling parameters of the model are: (1) the 
local stripping rate of the liquid jet, M , (2) the local mean drop size 
being produced when is ‘•tripped from the jet, D*, (3) the droplet 
heating and burning rates, (4) the droplet drag coefficient, and (5) for 
the chamber flow, the rate of mixing of the external "Rigimesh flow" (for 
low percentage flow, this parameter will be shown to be of little importance). 

A correlation for the droplet drag coefficient was utilized which was 
obtained for accelerating, burning droplets in a convective flow field 
(Ref. 5). While it lacks particle interaction effects, it appears to be 
fairly accurate. 

The rate of mixing of the "Rigimesh flow" is important only if that flow 
is abnormally large. This mixing of the "Rigimesh flow" is not accomplished 
by turbulent mixing of striated parallel gas flows but is primarily 
"caught" or entrained between adjacent elements. As combustion of the spray 


proceeds, the reacting flow field expands radially as it progresses 
downstream as computed from the entire interrelated set of conserva- 
tion equations describing each flow field. The assumption of uniform 
pressure at each axial location is utilized to iteratively solve the 
equations. These iterations, coupled with the proportioned chamber area 
constraint, determine the pressure level (pressure varies axially but 
not radially), velocity and area of each flow field. Computations have 
been performed for extremes in "Rigimesh flow" mixing from complete mix- 
ing to no mixing. Results show that rapid spreading of the combusting 
flow field from adjacent elements (with the presence of "Rigimesh flow", 
mixing or not) causes a decreasing axial pressure. The "Rigimesh flow" 
accelerates rapidly so that, within approximately 2 inches downstream from 
the injector face, its flow area is reduced to only approximately 3 percent 
of its injection area. Because this area is in the form of an annulus 
trapped between elements, the average thickness of this annulus (for typical 
injectors) is on the order of 0.01 inch. Thus, the "Rigimesh flow" allows 
a single turbulent eddy to sweep the flow into the adjacent element flow 
fields. Consequently, with normal amounts of "Rigimesh flow" (approxi- 
mately 5 percent of the total fuel flow), the axial variation of the 
expansion area of the combusting flow field of adjacent elements is 
relatively unaffected by the presence of the "Rigimesh flow". Hence, r.he 
rate of mixing of the "Rigimesh flow" is not of importance and is usually 
taken to be a linear function between the injector face and the 95-percent 
closure point of two adjacent elements computed under no Rigimesh mixing 
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conditions. If this approach requires too much computer time, an 
arbitrary downstream distance may be chosen with no significant change in 
accuracy. Naturally, when the axial rate of mixing is prescribed, the 
mixed "Rigimesh flow" is spread uniformly over the cross-sectional area 
of the element's flow field and becomes part of the fuel to be reacted. 

An important aspect of this program was the development of a continuous 
sub- and supercritical burning (and heating) rate drop model that allows 
steady-state combustion analyses and performance predictions to be made 
to >5000 psia. The equations are similar to the El Wakil equations (Ref. 6). 
However, the boundary condition was changed to allow the existence of an 
external mass flux (i.e., surface regression effects of the droplet). 

This change allows smooth computation through the critical point. 

Additionally, high-pressure effects due to the presence of other gases 
were included for the computation of the vapor surface mole fraction and the 
"heat of vaporization". This involved use of the Redlich-Kwong equation 
of state, fugacity relationships, and solubility effects of the external 
gas in the droplet (Ref. 7) . 

Equations for the jet stripping rate and drop size production are pre- 
sented below: 


A. Stripping Rate 


V C A 


M j (P * U r )2 

'FT 


11/3 


TTD^AZ) 


Atomization 

Coefficient 


B. Mean Drop Size 


D. 

J 



U i (g / p 3 ) 




Drop Size 
Coefficient 


1/2 


2/3 


where 


Dj = diameter of jet 

U = expanded relative gas velocity in cup (or chamber) 
r between gas and spray 

Z = axial location 

p = gas density 

= jet density 

u. = jet viscosity 
J 

o. = jet surface tension 


The equations include an atomization coefficient, C^, and a drop size 
coefficient, B . It is not expected that the values for and C A be 

i\ 

the same inside the cup region as out in the chamber since the fuel is 

not constrained in the chamber as it is in the cup. Thus, the values of 

B and C in the cup may reflect droplet breakup, etc. However, the 
A A 

equations are general in nature and the values of and C A determined 
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for either regime have been found to be the same for any engine operating 
condition. This is significant in that the extremes range from burning 
to non-burning propellant cup conditions, various chamber shapes, etc,, 
and even different propellants. 

Methods of evaluating and for the cup and chamber regions are 
required and the validity of these equations has been verified. 

and C A for the cup region were determined by analytical comparison of 
four to seven different cases where the pressure drop from the end of the 
fuel gap annulus to the injector face (cup pressure drop) had been measured. 
Four of these cases were from subscale firings of LOX/hydrogen coaxial 
injection engines: (1) the SSME straight oxidizer post preburner, P c = 

1500 psia; (2) the SSME tapered oxidizer post preburner, P £ = 1500 psia; 

(3) the stability preburner- like uni-element motor, P c = 500-1000 psia, in 
which the delta P between the fuel gap annulus static pressure and the 
injector face pressure was measured directly; and (4) the SSME uni-element 
tests, P c = 1500 psia which consisted of a uni-element preburner and a uni- 
element main injector containing an oxidizer post capable of being recessed 
from flush mounting to 0.3+ inch from the injector face. The remaining 
cases were full-scale firings of three similar engines also using L0X/GH 2 
propellants: (1) the J-2, (2) the J-2S, and (3) a variable-oxidizer-post- 

recessed aerospike (segment) engine, P £ ■ 750 psia. The J-2 and J-2S are 
the only ones in which the propellants ignite and burn within a flared 
fuel cup sleeve. 
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The number of experimental cases used to determine B A and C A depended 
on the method used to compute the pressure profile from the fuel gap 
annulus to the point of propellant contact. At least two sets of data 
are needed to determine and C^; these resultant values are then used 
to predict the cup APs of the remaining engines. The adequacy of the 
chosen B a and C A is determined by comparing predicted and measured cup 
AP values. Acceptable accuracy has been taken to be ±20 percent (except 
at very low cup AP's where predicted differences of only a few psi can 
result in apparent high percentage deviations) . 

Three methods have been used for computing total cup delta pressures 
and for obtaining B and C.. This computation concerns the process 

i\ /\ 

through which the gas is assumed to flow from the gap annulus into the 
recessed portion of the injector element. In the first of these methods, 
the gaseous fuel was assumed to flow around the pest thickness and fill 
the entire annulus between the liquid jet and the fuel sleeve walls. The 
liquid jet was not allowed to expand, and therefore does not have the 
bulge as shown in Fig. 5. The fuel pressure in the annulus gap is calculated 
by assuming that the fuel static pressure did not change in the expansion 
from the annulus gap to the cup. B A and C A for the cup region were 
determined by comparison of predicted and measured cup AP values for the 
J-2 and J-2S. 


Sets of values of B and C were calculated that yielded the correct 

A A 

pressure drop for each engine. When these were plotted, they produced a 

cross point of two curves that was sharp and yielded the only values of B A 

and C that satisfied the pressure drop of both engines. It was necessary 
A 

for these calculations to assume burning to occur within the cup. A 
nonburning assumption caused over 90 percent of the jet to be atomized in 
the 0.200-inch recess of this cup to yield the correct pressure drop, 
clearly this was a nonrealistic assumption. Additional engine cup AP data 
were then used to verify the validity of the chosen B^ and C A> Although 
this method of obtaining B A and C A (combined with the fact that ignition 
of the flowing gas was assumed to occur at a cross-sectional plane just 
downstream of the oxidizer post) gave satisfactory agreement with most 
existing cup pressure drop data, it does not properly conserve gas momentum 
and, therefore, was considered unsatisfactory. Further, this technique 
failed to properly predict the directly measured cup AP's of the stability 
preburner like-element injector. 

The second method involved the use of an iterative procedure to predict 
a "sudden expansion" of both the gas and liquid (requiring the equations 
of motion for both propellants to be satisfied) with compressibility 
factors being used to allow for nonideal gas behavior. The iteration was 
used to find the allowable areas of expansion for the two propellants 
subject to the constraint of the total cup flow area. As expected, the 
results indicated that as the thickness of the LOX post is increased the 
pressure recovery approaches zero. Again values for B A and C A were 
obtained from the J-2 and J-2S data allowing ignition to occur as in the 


previous case. 


It became evident at this point that ignition of the cup gases upstream 
of the flare was not possible because of the high gas velocity in the 
unflared section of the cup. However, assuming ignition to occur at a 
cross-sectional plane located at the beginning of flare, failed to 
produce a common solution when the B A 's and C^'s were plotted. (The 
assumption that ignition occurred in the flare was also used to obtain a 
set of values for B A and C A with the first method but the change in the 
results was insignificant.) Although the proper values for B A and C A 
could not be adequately defined by this second method, comparisons were 
made to other engines. Using this approach the measured subscale SSME 
preburner total manifold-to-injector end pressure drop was accurately 
predicted. However, when predicted pressure drops were compared with 
other available large engine cup AP data the comparison was again less 
than satisfactory. Further, the trends of the directly measured cup 
pressure drop of the subscale preburner- like unielement firings could 
not be properly predicted with this method. 

The third method corresponds to a significantly different concept of 
the flow behavior within a recessed oxidizer post coaxial injector. The 
previous methods considered the gas and liquid static pressures to be 
equal (and radially uniform) at the end of the oxidizer post. Consequently 
the radial liquid and gas expansions predicted by the second method 
always occurred within the radial thickness of the oxidizer post. 

However, base bleed analyses of the injection process with, e.g., a 
coaxial element, suggested that the fuel pressure and oxidizer pressure 
(and the intermediate pressure at the tip of the post) may not be equal 
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at the end of the post and that adjustment of the streams to reach a 
radially uniform pressure would occur. Further investigation led to the 
concept previously illustrated in Fig. 5. The liquid jet pressure at 
the end of the oxidizer post is believed to be initially less than that 
of the surrounding gas. Relaxation of the radial pressure profile to 
a uniform pressure occurs due to initial expansion of the liquid jet 
flow area and subsequent contraction of the fuel flow area. Consequently 
the liquid jet flow gains static pressure and the gas flow loses static 
pressure until the pressure is radially uniform. In the analysis this 
initial liquid expansion is allowed to occur in approximately one 
oxidizer post thickness downstream from the post. Propellant contact 
begins at the point marked 1 in Fig. S and subsequent atomization of the 
liquid begins. Total cup AP is still defined as the difference between 
pressure in the fuel gap annulus and injector face chamber pressure. 

To compute I$ A and C A for this third method at least four sets of engine 
data are needed, because the equations used to predict the liquid area 
expansion also contain two empirical constants. Actually, data from 
all of the engines were used to obtain a "best-fit" correlation for 

injector cup pressure losses under non-burning cup conditions, as shown 
on Fig. 6. 
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Figure 6. 


Injector Cup Pressure Loss Correlation 
(Nonburn I ng Cups) 




Conservation and correlation equations were used to calculate the 
liquid jet area (and the contracted area of the fuel). Because the 
manifold-to-injector-face pressure loss for the fuel has the most 
significant effect on engine design, the conservation equations (mass, 
momentum, energy, and state) for the fuel were applied to calculate 
this fuel flow field from the end of the post (i.e., from the fuel gap 
annulus) to the point of propellant contact. The area of the contracted 
gas flow was correlated the following equation: 


A 

— EE. 

A ANN 


0.982 ♦ 0.0337 Ln f , p0St I,P ‘ 


g 


fuel sleeve 


where A is the area of the gas at the initial stream contact point 

and A^^. is the area of the fuel gap annulus, 
g 

In addition, the constraint was used that the sum of the areas of the 
gas and liquid at the contact point equaled the total area of the cup 
(fuel sleeve). Simultaneous with the determination of the two constants 
in the above equation, the cup and values were also determined by 
computing cup pressure losses downstream from the propellant contact 
point. 


The B a and C A values for this were determined from non-burning cup 
gas data. Resulting were the first values for cup B A and C A that yielded 
predicted pressure losses in agreement with data from the stability uni- 
element injector. Because the values for B A and C A were obtained 
assuming a non-burning cup, the J-2 and J-2S cup pressure losses were 
computed by assuming cup gas ignition to occur at a conical surface 
based at the upstream edge of the flare on the fuel sleeve and sloping 
downstream toward the liquid jet. Also, because B A and C A and the J-2 
and J-2S pressure losses are all known, this method allowed the flame 
speed in this reacting two-phase LOX/hydrogen flow to be predicted. The 
value obtained for the J-2 was 382 ft/sec. The fiame speed for the 
J-2S was predicted to be slightly greater than its gas velocity at the 
beginning of the flare. Consequently, the latter flame speed was reduced 
to this gas velocity. As a result, the predicted cup pressure loss 
for the J-2S aoes not precisely match the measured loss (as shown in 
Fig. 6). Nevertheless, this third method is preferred and is recommended 
over the first two methods. It is believed that as more accurate 
non-bunting cup gas pressure loss data becomes available, it will be 
possible to adjust B A and C A to allow prediction of the measured loss 
within ±10 percent. 
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Although a nonplanar ignition front has been introduced, the program is 
still one-dimensional. Flame propagation is computed using an 
"averaged" or pseudo mixture ratio in the gas flow field region 
surrounding the liquid jet. This region is composed of unburned 
propellant (gas flow near the jet that has not yet passed through the 
"flame front") and burned propellant which has passed through the flame 
front. To compute the mixture ratio under such quasi-two-dimensional 
conditions the following equations were utilized. 


y l ' y f(x) 


I l 2 2 

\ y I ' y jet(x) 


where 


r f(x**x) * y f(x) " Vj 


x * axial distance 


MR, . - the useudo mixture ratio at axial location x 
(x) 

MRj > the mixture ratio at the ignition point 
/j ■ the radial distance from the centerline of the liquid 
jet to the ignition point on the fuel sleeve 
y^ x ^ ■ the radial distance from the centerline of the liquid jet 
to the flame front at location (x) 


Xi ./ * ■ the radial distance from the centerline of the liquid jet 
7 jet(x) 

to the outside surface of the jet at location (x) 
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V^. ■ the flame speed 

Vj ■ the gas velocity at the ignition point. 

Should part of the gas flow at the injection point be composed of 
fuel rich combustion products (i.e., topping cycle engines) then 
appropriate adjustment to the mixture ratio equation is performed by 
the program. The program uses the local mixture ratio (MR. ) with 
the combustion gas properties table to calculate corresponding pseudo 
gas stagnation temperature, etc., to the local pressure through interpo- 
lation techniques. This stagnation temperature is then reduced in the 
program by the corresponding amount of energy stored within the 
unreacted oxidizer vapor and droplets at location x. Static temperature 
is computed assuming "frozen*' composition as will be described later. 

Values for and are also needed for the chamber. To obtain these, 

comparisons were made with data from two different engines. In this 

case the comparisons were made of measured and predicted combustion 

efficiencies of short chamber length, non-burning-cup segment engines. 

Utilizing the third method developed for the cup, curves of B and C 

A A 

which give the correct performance level for each engine were obtained. 
These were obtained by calculating the gas and spray flow and combustion 
in the chamber. Nonbuming cup exit flows were ignited ?t the cup 
exit on the fuel sleeve and a flame speed of 600 fect/sec was utilized. 
Flame speed in the chamber is not precisely known, however, parametric 
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studies show that flame propagation is complete in no more than .15 
inches (three normal axial steps) for flame speeds as low as 350 feet/ 
second. Cross plots of B and C for the proper chamber conditions 

ft 

yielded a single cross point. The final point selected was checked 
by predicting an independent segment engine efficiency and also 
predicting the liquid jet lengths of the 3-2 and J-2S. Combustion 
within the latter two engines had been observed and photographed in 
subscale transparent hot-firing tests (Ref. 8). The model predictions 
are in good agreement with the test results when the computations are 
performed without using DER/STC but using the CICM program to compute 
combustion efficiency q # down to the engine throat. Computation 
of B. and C for the chamber was performed in a similar manner since 

r\ 

the segment engines used to obtain the values of the parameters were 
especially chosen. None of the engines had Rigimesh flow, all had 
100 percent efficient intraelement mixing characteristics, and no 
manifold flow maldistribution. Consequently the engines could be 
viewed as single streamtubes in which the mixture ratio was uniform. 
Additionally a great deal of measured performance efficiency data 
were available for these selected segment engines. 

SUMMARIZATION OF THE CONSERVATION EQUATIONS FOR SINGLE ELEMENT 
COAXIAL INJECTION COMBUSTION 

As stated previously the program is based on the use of conservation 
equations for both the liquid jet and combustion gas/spray flow in the 
cup region and chamber flows. Where appropriate, as in the cup, these 
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equations have been used with other correlation equations to determine 
the propellant contact point and the flow conditions of the propellants 
at that point. In addition, the effects of flame propagation in flared 
cups and in the chamber (from flow exiting from non-burning cups) 
have been introduced. Values for and C A in the cup and chamber have 
been determined. 

The analysis of the injector element, as far as prediction of vaporiza- 
tion efficiency is concerned, proceeds from a flow field such as that 
depicted in Fig. 5. Three axisymmetric flow fields must be considered; 
(1) the liquid jet, (2) the spray/gas burning (or non-burning) flow 
surrounding the jet; and (3) the Rigimesh flow (in the chamber) 
surrounding the spray/gas flow and mixing with it. Flow within each 
of these fields is considered to be one-dimensional. 


Four sets of conservation equations art used, one each for the 
combustion (or non-burning) gas, the spray, the liquid jet, and the 
Rigimesh flow. The four sets of equations are related by expressions 
which describe the transport phenomena between the flow streams. 

It is quite possible and convenient to sum the conservation equations 
for each constituent and obtain one overall set of combined equations 
replacing one of the four original sets of equations. Numerical techniques 
and practice indicate that this set of combined equations should replace 
the combustion gas conservation equations. Terms for transport phenomena 
between streams are decoupled from the combined conservation equations 


so that they can be directly integrated. The entire set of all 


equations (and their derivations) can be found in Ref. 8. For this 
user's manual only a verbal description of the conservation equations 
are presented. However, the expressions representing the transport 
phenomena of mass, energy and drag force are presented in detail. 



A. 


The Overall Combined Equations 

i 

1. The local mixture ratio equation. 

2. The overall continuity equation expressing conservation 
of the mass of the gas, spray, jet and Rigimesh flow at 
every increment. The sum of the mass of each flow is 
a constant . 

3. The overall momentum equation 

4. The overall energy equation 
a) This equation can be rigorously written for the flow 

enthalpies and velocities and with the use of 
tabulated gas equilibrium properties as a function of 
mixture ratio and pressure solved iteratively with the 
other equations. This requires extensive triple 
interpolation and produces an inefficient program. 

Experience indicates that the equilibrium combustion 
gas properties are weak functions of the stagnation 
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pressure; hence the overall energy equation is 
replaced by a set of combustion stagnation properties 
dependent only on the injection pressure and tempera- 
ture and the local (axial varying) mixture ratio. 
Stagnation temperature, etc., (as a function of axial 
location) is then directly and easily computed from 
the properties tables. This temperature is adjusted 
(reduced) for the energy contained in the spray and 
remaining oxidizer vapor. Reduction to static flow 
temperature is accomplished by assuming that no 
species change takes place between stagnation and 
flow conditions. Thus, except for temperature, gas 
properties (specific heat ratio, viscosity, etc.) 
in the flowing stream are assumed to be the same as 
at stagnation. The equation 
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is used to compute the static temperature, 
Here T ■ static temperature 

T q * stagnation temperature 

Y Q ■ specific heat ratio 

MW q ■ stagnation molecular weight 
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R = universal gas constant 
V = velocity of gas 


Good agreement results when this method is compared 
with the rigorous equation (and its attendant 
compacted properties tables). 

5. Equation of state 

B. The Jet Equations 

Expressions for continuity, momentum and energy are considered; 
drag on the jet is neglected in the momentum equations 
since the effect is accounted for in the production and 
acceleration of droplets. Similarly, the jet temperature 
is considered to be constant since the surface stripping 
prevents conduction to the jet core. 

C. The Rigimesh Flow 

This flow is considered to be isentropic. The continuity 
equation contains a mixing rate expression, but this does 
not affect the conditions required to have isentropic flow. 

D. The Spray Equations 

It is this set of conservation equations, along with the 
expressions for the stripping rate and dropsize production 
tha. essentially control the program. Of particular importance 


in the spray equations is the vaporization rate, drag force 
and heating rate expressions in the continuity, momentum and 
energy equations, respectively. These equations are principally 
initial value problems in that a new initial condition is formed 
at each axial increment along the jet. 

All of these conservation equations are required for simultaneous 
solution on a digital computer to predict engine performance. Iteration 
of the initial assumed injector face pressure is required until the 
throat velocity is sonic. The computation of the combustion gas is 
considered to be composed of constituents in thermodynamic equilibrium. 
This is in agreement with the accepted approximation that, for well 
designed engines, drop vapor diffusion rates are very much more 
limiting than gas phase chemical kinetic rates. 


Expressions Describing the Transport Phenomena 
1. Drag Force on Droplets 

The expression describing this droplet dynamic transport term 
appears in the spray momentum equation. The drag force is de- 
fined as: 


F ° " 8g 


p (D n ) 2 (u - u n ) ( | u - u n | ) C n 
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- droplet drag coefficient, 
initial drop group size n. 


n 

I) 

f ,n 


P 

u 


X 

p 


- droplet diameter 

- drag force on droplet 

- gas pressure 
-gas velocity 

- droplet velocity 

- axial location 

- gas free-stream density 


The drag force includes both frictional drag and the drag due to volume 
forces across the drop arising from any existing gas pressure gradients. 
Other terms in the drag force equation, such as the acceleration of the 
"apparent mass" of the gas displaced by the droplet and the Basset term 
(non-steady condition) have been neglected because 

p n >> p of the gas 

p n - the density of the droplet 


The validity of the equation is limited to the applicability of the existing 
correlations for the droplet drag coefficient. The review of existing 
correlations in Ref. 8 indicate that Rabin's, et al, work in Ref. 5 is 
still considered to be the best correlation for describing drag coefficients 

s' 

| when ;» pp I i pd to drop I ft « in a rocket combustion chamber flow. Rabin's work 
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shows that c^' 1 is a function of the relative Reynolds number. The 
correlation includes (1) the effect of gassification in a convective 
flow field and the effects of distortion of the drop. 


C D = 24 Re' 0-84 Re <. 80 

= .271 Re 0,217 80 < Re < 10 4 

=2 Re > 10 4 

where 



Re = 

g v 

and u is the gas free-stream viscosity. 

2. Droplet Vaporization and Heating Rate* 

Background . The quasi-steady evaporation coefficient approach 
to droplet heating and burning, while empirically based on the 
observation that a burning droplet's diameter squared varies 
linearly with time, has been expressed analytically in increas- 
ingly comprehensive formulations. These models are based on the 
concept that a spherical flame surface surrounds a spherical 
droplet, with simultaneous heat transfer to and evaporation from 
the droplet being enhanced by the presence of the flame. These 
models have all been formulated as quasi-steady problems (i.e., 
time variation has been neglected in writing the conservation 
equations) , although there are no assumptions in the models 
that preclude droplet heating. Relatively recent work at 

*For clarity, notation has been changed in this section; nomen- 
clature is placed at end of section. 
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Rocketdyne (using the addition of diffusion equations) has culminated 
in the added development of a thin- flame model that includes uniform 
droplet heating. A problem that arises in applying such a model, how- 
ever, is that the initial heating and burning rates may be over-pre- 
dicted by assuming a flame exists when the vapor concentrations are 
too low to support it. Another problem is that the derived formulae 
for the burning rate (or the evaporation coefficient), in all of these 
models have singularities (blow-up logarithmically) if droplet tempera- 
tures approach propellant critical temperatures. One final problem is 
that exposing the droplet to even mild forced convection is likely to 
blow the flame into the droplet wake or extinguish it, so that flame- 
enhancement of vaporization does not occur. 

As a consequence of these limitations and problems, propellant droplet 
gasification and burning has also been analyzed from a vaporization 
standpoint, with vapors diffusing into and mixing with a high-tempera- 
ture gas stream. So far as the droplet is concerned, combustion reac- 
tions within that gas stream serve to keep the gas temperature high and 
the vapor concentration low. (In practice, reaction to local thermo- 
dynamic equilibrium is usually assumed.) To the extent that the free- 
stream gas temperature is lower than the stoichiometric flame tempera- 
ture (the thin-flame model's driving temperature), a vaporization model 
will predict lower droplet burning rates than will a thin-flame model. 

An evaporation model that is commonly used for analyzing spray gasifi- 
cation in rockets is that of El K».kil (Ref. 6) and others at the 


University of Wisconsin. By solving spherically symmetric, quasi- 
steady conservation equations for simultaneous heat and mass transfer, 
droplet mass evaporation rate and (uniform) heating rate expressions 
have been developed. 


It is possible to calculate nonuniform temperature distributions within 
a droplet undergoing heating (e.g., Ref. 8), but it is usually assumed 
that internal temperature gradients are prevented from building up by 
strong internal circulation. Under convective flow conditions, surface 
shear does promote circulation and this simplification is probably quite 
valid. Then the uniform droplet temperature is obtained from: 


d(T n '~l 

r 3 1 

r n . n n, 

" r 

6 P £ ( p£ V D ) 
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Forced convection and resultant nonspherical transfer processes are 
accounted for through empirical Nusselt number correlations for both 
heat and mass transfer. The Nusselt number correlations used in the 
mass transport equation were obtained by Ranz and Marshall (Ref. 9); 
based on droplet film (f) conditions. 


Nu » 2 (l + 0.3 Sc. 1/3 
m v r 

- 2 (l +0.3 Pr f 1/3 


* f 1/2 ) 

< /2 ) 


They verified this correlation with data from vaporization of water 
droplets in heated air. The equations derived thus account for both 
droplet heating and evaporation. 
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The foregoing droplet heating and evaporation model is capable of 
computing droplet behavior to complete combustion at subcritical 
chamber pressures, although the vaporization rate blows^ip loga- 
rithmically as droplet temperatures approach the boiling temperature 

(X‘' -► 1). For most conditions, the "wet bulb" effect suppresses the 

v d 

equilibrium droplet temperature enough below the boiling point to 
avoid the singularity. There, however, the evapoiation rate is 
strongly dependent upon droplet temperature and, because an implicit 
solution of the system of equations is required, many iterations may 
be needed to obtain convergence. Recent work, summarized in Ref. 8, 
gives good correlation with experimental data under such conditions, 
even up to high pressures, if the effects of the presence of other 
gases on the vapor pressure and "heat of vaporization" are taken into 
account . 

Real Gas Effects . For vapor-liquid equilibrium, the free energy is the 
same on either side of a phase interface. This fundamental relationship 
for vapor-liquid equilibrium is conveniently written in terms of 
fugacities; for each component i, the fugacity of the vapor f/ is equal 
to that of the liquid f^ (Ref. 7). Since the liquid senses the total 
pressure while the vapor senses only its partial pressure, the equilibrium 
relationship can be written as 


f i ( P v) " f i ( P Total) 




Hence, at constant temperature, as the total pressure increases, the partial 

pressure of the vapor has to increase to maintain the required relationship 

for equilibrium. For a non-ideal gas, the entha.py is a function of the 

partial pressure of the gas (Ref. 10). Hence, the heat of vaporization, 

AH will be a function of total pressure since 
vap’ 


In the calculation of vapor-liquid equilibrium, the vapor has to be considered 
a non-ideal gas. Of the four two-constant equations of state which have been 
widely used, the Redlich and Kwong equation is more accurate and the best at 
high pressures. The Redlich-Kwong equation is. 

RT a 

= Cv-b5' " t °-5 v(v+b) 

where a and b are determined from mixing rules (Ref. 7). To match data over 
wide ranges, a and b have been programmed as functions of temperature. 

These "real gas" corrections have been neglected in most prior applications of 
the El Wakil droplet vaporization model. Under supercritical pressures, some 
conditions led to calculated equilibrium temperatures below the critical 
temperature, but usually no equilibrium temperature was reached and the drop- 
lets were heated through the critical temperature. The model could be used 
beyond this point, but it usually was not because a physical model was lacking 
for xj at the "surface" of the pure supercritical vapor pocket. Instead, most 
users either assumed instantaneous mixing of such supercritical vapors with 
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the surrounding gases, which is obviously unsatisfactory, or switched to a 
supercritical burning model due to Spalding (Ref. 11). This latter model, 
however, treats only the mass transfer and assumes that the vapor pocket 
remains at its critical temperature. As a result, no prior combustion model 
employin the El Wakil vaporization formulation can be adopted carte blanche 
for supercritical spray heating and combustion. 

Interestingly, introduction of the real gas corrections for vapor pressure 
and heat of vaporization caused the El Wakil solution for droplet temperature 
to reach a subcritical equilibrium temperature for all conditions. This is 
known from photographic evidence (Ref. 12) to be unreal, so the need for an 
improved formulation was apparent. 

New Droplet Heating an d Diffusion Model . The El Wakil model has been extended 
and improved to overcome this physically unrealistic result. The new model 
is referred to as the droplet diffusion model. The main difference between 
it and the old model is this: In the El Wakil formulation, only the propel- 

lant vapor is considered to have a non-zero net flux in the film surrounding 
the droplet, while in the new model the radial mass flux of combustion gas in 
the film surrounding the droplet is no longer assumed to be equal to zero. 
Instead, the molar flux of combustion gas is defined at the droplet surface 
through a moving control volume formulation such that changes in the droplet 
radius, due to droplet density changes and mass diffusion, cause it to be 
greater than or less than zero. The droplet surface boundary condition is 
determined through use of the species continuity equation. This is one of 
the major changes developed since the initial version of this model was pro- 
grammed into the droplet heating and vaporization version of the current 
DER(STC) computer program. (The other major change is the inclusion of 
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solubility, us ins the methods of Ref. 7, of the external gas in the droplet 
surface layers; this latter change allows computation of the droplet surface 
vapor concentration to extremely high pressures.) The droplet surface boundary- 
condition equation? are: 


a) 

MW_ N_ = 

n 3r" 

p_ d 

E d E d 
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Thus, as the droplet "burns' the external diffusing combustion gas is allowed 
to enter the control volume and occupy that fraction of the volume vacated 
by the receding droplet surface. 


The diffusion rate, or burning rate, is defined by the diffusion equation and is 



where 


and 


B = 
(NOTE: 



Here f refers to "film" conditions.) 
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The droplet heatup rate is defined to be 


Q n = TTk f Nu h D n z 

where 

•n 
m 

v 

TTk f D n NUj^ 

The droplet diffusion model no longer has the logarithmic singularity at 

either the droplet boilinp or propellant critical temperatures because, as 

droplets are heated through these temperatures, the value of B is such that 

(.1 - Bx 11 ) does not vanish. It thus becomes possible to continue analyzing 
v d 

spray droplets' behavior after they have become fully gasified, but have not 
yet been diffused and mixed into the surrounding combustion gas stream. 



i 

i 

\ 

! 




Comparison of the foregoing droplet diffusion model equations with the old 
model equations, e.g., as given by El Wakil, shows them to be very similar. 
The major differences are the appearance of the parameters A and B. Examina- 
tion of the equations shows, however, that A and B depend upon the heating 
and vaporization rates so that the droplet diffusion model must be solved 
implicitly by iterative methods. If the heating and vaporization rates are 
low enough that 3r”/3t vanishes, A+1 , B-*T and the droplet diffusion model 
reduces rigorously to the El Wakil model. Chemical reactions are not taken 
into account directly in the droplet heating and diffusion model, but combus- 
tion is simulated by specifying a bulk gas equilibrium flame temperature 
and zero droplet vapor mass fraction in the local free stream (except where 
a flame front is around the jet). 
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For this section: 
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parameters in Redlicn-Kwong state equation 

parameter in droplet diffusion model 

specific heat at constant pressure 

molecular diffusivity 

droplet diameter 

fugacity 

enthalpy 

heat of vaporization 

thermal conductivity 

molecular weight 

rate of change of mass 

the absolute gas molar flux 

Nusselt number 

pressure 

Prandtl number 

spray or droplet heating rate 
Reynolds number 
universal gas constant 
radial coordinate (drop radius) 

Schmidt number 

temperature 

time 

molar or specific volume 
mole fraction of droplet vapor 
heat transfer blockage term 
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Greek Letters 


p density 

Superscripts 

L liquid 

n concerned with tiie droplet size group 

v vapor 

Subscripts 

d droplet (droplet surface) 

E external gas 

f droplet film 

h heat or heating 

l liquid (usually referring to droplet properties) 

m mass 

v droplet vapor, vaporization rate 
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INTERFACING WITH DER 

Integration of the computer program into the DER methodology was formulated 
after considering many alternative methods for handling intra-element mix- 
ing. It was decided that the most accurate way to include intra-element mix- 
ing was to divide the spray and gas flows for each coaxial element zone 
into multiple mixture ratio zones. The manner in which the zones are sub- 
divi, . to simulate the intra-element loss is determined from cold flow 
measurements. These cold flow measurements relate the element geometry and 
flow conditions to its mixing efficiency. Such information is input to the 
CICM program in terms of mass fraction as a function of the total fuel and 
oxidizer flowrate fo. each coaxial element zone. An auxiliary program (sub- 
program) is used to perform this analysis and interface CICM with the stream- 
tube portion of DER. The auxiliary program organizes the spray/ gas informa- 
tion generated by the coaxial element zone calculations in terms of punched 
card output. The strcamtube portion (STC) of DER must, of course, be pro- 
vided with more information than this punched card output. 

During a coaxial element zone CICM analysis, as many as 100 droplet groups 
can be generated by the stripping process. However, to interface with the 
SIC portion of DER, the droplet groups must be condensed to fewer than 12 
equivalent droplet groups (restricted by STC program). The auxiliary program 
condenses the number of CICM droplet groups to those necessary through 

use of input variables which define the mass fraction of the spray in each 
DER droplet group. The DER droplet group temperature, velocity, and droplet 
diameter are determined by requiring conservation of droplet energy, droplet 
momentum, and droplet spray vaporization rate, i.e., 
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where the lower and upper limits of summation are determined by 


IV = f 
DER. SPRAY 

J DER 


AY f W SPRAY ) 

■ . \ TOTAL / 


The amount of fuel and oxidizer for each DER zone are specified as mass 
fractions of the total fuel and total oxidizer flowrates for each coaxial 
element zone calculation. The DER start plane pressure, zone areas, and 
gas velocities are determined by requiring conservation of the gas mass, 
momentum, and energy in the coaxial element zone with the constraint that 

the sum of the DER zone areas must be equal to the chamber cross-sectional 
area. 

The output of the auxiliary program consists of streamtube initialization 
punched cards for use in the supercritical version of the DER program. 


The order in which the st reamt ube initialization cards are punched depends 
on the order in which the coaxial element zone calculations were performed. 
The first set of streamtube cards corresponds to the first coaxial element 
zone calculation, the second set corresponds to the second coaxial element 
zone calculation, etc. The user can, if he chooses, reorder the streamtube 
cards in any manner that he selects before executing the DER program. 


SO 


MAIN PROGRAM 


A logic diagram of the CICM main computer program is shown in Fig. 7. 

Also shown in Fig. 8 through 11 are logic diagrams for several important 
subroutines. The method of solution used in the main program is summarized 

in the following paragraphs. 


Input data required are: (1) tables of propellant and combustion gas prep 

erties, (2) properties of the equilibrium combustion gas at stagnation 
conditions, (3) miscellaneous program control information, (4) case infor- 
mation data. The input data are printed as they are road by the program, 
which permits a full documentation of the computer run conditions. 


The input data are used in an initialization section to calculate a number 
Of program variables which include updating the stagnation equilibrium com- 
bustion gas properties (CGTBI2) , defining the cross-sectional area as a 
function of axial distance (AVAR or AVARR) , and velocities and properties 
(IMR) at the start plane. Initialized data are printed out before entering 

the main computational iteration loop. 


The main computational loop solves the model iterations at each axial posi- 

tion using sequential marching numerical methods. At each axial 

the liquid jet stripping rate and mean droplet size generated by the strip 


ping are calculated based upon the local combustion gas velocity and combus- 
tion gas properties (ATM; a new spray droplet size group is initialled 
from these local data. Droplet acceleration (DRAG) , heating, and vaporization 
(DIIVS) are then calculated for each droplet group present in the combustion 
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Figure 7. CICM Main Program Flow Chart (cont) 
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Figure 8, Subroutine CGAT Flow 
Chart (cont) 
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Figure 11. Subroutine FDTDX Flow Chart 


























gas. A portion of the "Rigimesh" gas is then mixed into the combustion gas 
and downstream gas velocities and properties are calculated based upon the 
total droplet vaporisation rate, amount of "Rigimesh" gas added, and the 
cross-sectional area (CGAT) . Initialization of parameters for the next step 
is then performed (INIW) and, at selected axial locations, complete gas and 
propellant spray group data are printed (0UTPUT) . 

Upon completion of the main iteration loop, if the case was a cup calculation, 
the program then checks to see if the cup exit pressure is within a tolerance 
(which is input) of the chamber pressure. If the cup exit pressure is outside 
the tolerance, the case is rerun with a new estimated cup delta pressure. If 
the cup exit pressure is within the tolerance and the case is a coupled cup/ 
chamber calculation, the calculated cup exit conditions are used as initial 
conditions for the chamber calculation along with chamber information (INCIIA) . 
Also, if the DER option w'as specified in the input, cup exit conditions are 
saved on a scratch tape (0UTCUP). If the case is not a coupled cup/chamber 
calculation, the program branches to a read location specified by the variable 
IREAD to begin a new case or terminates the calculations. 

If the DER option was specified, spray and gas data are saved (0UTDER) on a 
scratch tape unit. If all the DER zones have not been executed, the program 
branches to the case input statements to read in new zone data. If all the 
DER zones have been executed, the program checks to see if each chamber case 
was continued to the axial position required for DER punched card output. If 
any of the zone chamber cases was terminated before reaching the axial posi- 
tion required for DER punched card output, the program recalculates these 
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, 0 „ e chamber cases f INCUP, . Upon completion #f aU :one calcu]ations ^ 

program calculates the DER peached card information (DERINI) , plmc hes the DER 

cords, and lists the DER punched card output. The program then branches to a 

read location specified bp the variable IREAD to begin a new case or terminates 
the calculations. 


The version of CICM describe in 

aescribed ln 'hie report consists of a main or calling 

program together with 33 subroutines. A listing of the CirM m 

g or tne CICM program, together 

»ith its subprograms and function routines, is shown in Appendix A. 



SUBROUTINES 


AT0M This subroutine calculates the portion of the liquid jet that 

is atomized over one axial computational increment. A droplet 
spray group is calculated, and the initial weight flow rate 
and initial mean droplet diameter of the group are determined. 

AVAR In this subroutine, the cross-sectional area per injection element 

is calculated at each axial computational step for an axisymmetric 
combustion chamber. 

AVARP In this subroutine, the cross-sectional area per injection element 
is calculated at each axial computation step for a combustion 
chamber specified by a table of areas at specific axial distances. 

CGAT In this subroutine, the portion of the "Rigimesh" flow which is 

mixed with the combustion gas stream is calculated. Combustion 
gas properties are reevaluated as a function of local mixture 
ratio. Division of the constrained area between the liquid jet, 
combustion gas, and "Rigimesh" streams are solved iteratively. 

CGPR0P In this subroutine, the combustion gas stagnation temperature and 
properties are calculated from tabulated values and the local 
mixture ratio of burned propellants. Droplet and non-reacted 
liquid vapor energies are subtracted from the tabulated stagna- 
tion temperatures. 
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CGTBIN 


CGTBI2 


CPLF 


CUBIC 


DERINI 


DHVS 


DINTRP 


In this subroutine, the stagnation equilibrium combustion gas 
properties are read into program as a function of mixture 
ratio. 

Subroutine entry point in CGTBIN which adjusts properties read 
by CGTBIN for changes in propellant inlet energy. 

In this subroutine, the liquid specific heat is obtained by a 
a double interpolation of values tabulated as a function of 
pressure and temperature. 

In this subroutine, the real roots, and the number of them, are 
determined from the coefficients of a cubic equation. 

This subroutine is used, if the option is selected for DER 
output, to recall data from a scratch data set. Spray droplet 
groups and gas flows are setup and parameters punched out for 
initial values to each stream tube in DER. 

In this subroutine, the heating and vaporization of each spray 
drop group are calculated with the support of several subroutines. 

In this subroutine, a linear double interpolation is performed 
using points and slopes which have been previously determined 
in subroutine L0CFAC. 
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DRAG 


EQSTAT 


FDTDX 


FGPR0P 


FT 


HEAD 


In this subroutine, the spray droplet velocity, which changes 
due to drag forces, is calculated. Droplet Reynolds number 
ane drag coefficient are calculated in the procedure, and 
constraints for the droplet velocity to approach gas velocity 
are imposed. 


In this subroutine, the densities of liquid and of gas mixture 
are calculated from the Redlich-Kwong equation of state. 

In this function subprogram, dT d /dx, is calculated for another 
function subprogram, FT. In the solution of d^/dx, spray 
vaporization and droplet diameter are also calculated. This 
function subprogram is basically the vaporization model in CICM. 

In this subroutine, mean droplet film physical properties and 
film diffusivity are calculated for subprograms FDTDX and DHVS. 

This function subprogram is used by subroutine DHVS, which also 
transmits it to subroutine ZER0, for determining the value 
of a function based on a predicted value of droplet temperature. 
This function goes to zero when the correct droplet temperature 
is found. The bulk of the calculations performed by FT are done 
in its subprogram FDTDX. 

This subroutine prints a header page to identify the computer 
program. 
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INCHA This subroutine establishes the chamber initial conditions 
following an injector cup analysis. 

INCUP In this subroutine, cup exit conditions, saved by subroutine 
0UTCUP, are read from a scratch data set. 

INDER In this subroutine, DER data, saved by subroutine 0UTDER, are 
read from a scratch data set. 

INIR In this subroutine, initial values for the flow parameters 

are calculated and printed. 

IN I W In this subroutine, parameters are initialized for next axial 

step. The "1" level parameters are set equal to the "2" level 
parameters of the previous step. 

L0CFAC In this subroutine, the location of the first of two sequential 
values in an array which bracket a specified value are found, 
and a scale factor, (X - X^/CXg - Xj), is calculated. The 
values in the array must be arranged in either ascending or 
descending order, and the validity of the order is checked if 
the option is specified. 

0UTCUP This subroutine causes cup exit conditions from an injector 
cup element analysis to be stored on a scratch data set. 


iZUTDER 


Output 


R1I0GF 


STLF 


TAB IN 


VISLF 


If the option is selected for DER punched output, this 
subroutine causes spray and gas flow data to be saved on 
a scratch data set for processing at the end of the job. 

This subroutine causes the solution at specified axial 
locations to be printed. This is the primary output routine 
of CICM. 

In this subroutine, the gas density is calculated. A com- 
pressibility factor is used, which is obtained by a double 
interpolation of values tabulated as functions of temperature 
and pressure. 

In this subroutine, the liquid surface tension is determined 
by performing a double interpolation of tabulated values as a 
function of temperature and pressure. 

This subroutine causes the propellant liquid and vapor physical 
property tables to be read into the program and printed out, if 
the option is specified. 

In this subroutine, liquid viscosity is obtained by a double 
interpolation of values tabulated as functions of pressure 
and temperature. 
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XVDHV 


ZER0 


In this subroutine, the mole fraction of vapor at the droplet 
surface, heat of vaporization, and Redlich-Kwong A and B 
parameters are obtained by interpolation of tabulated values. 

In this subroutine, droplet temperature is found, starting with 
upper and lower limits, by successive solutions of a function 
of droplet temperature with subprogram FT being used. The 
final value of the function must approach zero. If the primary 
numerical solution is unable to converge on a solution, a 
secondary numerical solution is automatically reverted to 
which has better numerical stability, but is less accurate than 
the primary numerical method. 


I’ROCRAM INPUT 


Specific input data fcr the CICM computer program are listed in 
Tables 1 , 2 , 3 , and 4 , which have been structured in the 
format of a typical input punched-card data deck. The input con- 
sists of blocks of cards describing the propellant and combustion 
gas, stagnation equilibrium combustion gas, control data, and case 
data. In these tables, the "CARD NO." is a suggested card identi- 
fication number (punched in columns 73-80) which is consistent with 
sequence numbers on the sample data cards listed in Appendix B. 

Where ranges of ID sequence numbers are given, consecutive integers 
are implied. (Note that the different blocks of the CICM program 
input data deck so sequenced should not be sorted with each other, as 

there is overlap and/or duplication of sequence numbers between these 
blocks . ) 


The "FORMAT" in the tables of input instructions denotes the type of 
FORTRAN input (integer, floating point decimal, alpha-numeric) and 
the subdivision of each card's first 72 columns into fields. Standard 
FORTRAN input formats are used. Specifically used are: 


Comment cards (A- format) 

Integer variables with variable names beginning with 
letters I through M (no decimal points, 12 space field 
widths, last digit in last space of field, 6 consecutive 
values per card until READ statement is finished). 
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TABLE 1 . INSTRUCTIONS FOR PROPELLANT AND 
COMBUSTION GAS INPUT DATA 


a :yj xc 


10 

(6112) 


VARIABLE 

CODE 


IPTAB 


20 

(6112) 


30, 31, etc. 
(6E12.8) 

40, 41, etc. 


50, 51, etc. 
(6E12.8) 


150, 151, etc 
(6E12.8) 


250, 251 , etc 


260, 261, etc 
(6E12.8) 


570 

(6112) 


NPTP(l) 

NPTP (2) 
NPTT(l) 

NPTT (2) 

TP ( I , 1 ) 

1=1 ,NPTP(1) 

TT(I.l) 

1=1 ,NPTT(1) 

(TXV (I , K, 1) 
1=1 ,NPTP(1) ) , 
K=1,NPTT(1) 


(TDHVCI.K.l) 
1*1 , NPTP (1)) , 
K=1,NPTT(1) 

TA(I,1) 

1=1 , NPTT (1) 

TB (1,1) 

1=1 ,NPTT(1) 


580,581 , etc 
(6E12.8) I 


NPCP (1) 

NPCP(2) 

NTCP(l) 

NTCP(2) 

TPCPL(K.l) 

K»1,NPCP(1) 


DESCRIPTION 


Print control integer: "0" to suppress printout, "1" 

to print out table of propellant properties.* 

Cards 20-264: input oxidizer droplet mole frac- 

tion, heat of vaporization, ana Redlich-Kwong 
equation of ctate parameters. 


Number of pressures in oxidizer mole fraction table. 

Limit: 2 to 30. 

Required to be zero(O) 

Number of temperatures in oxidizer mole fraction table. 

Limit: 2 to 20. 

Required to be zero (0) 

Pressure array in ascending order for oxidizer table. 
Units: psia 

Temperature array in ascending order for oxidizer table 
Units: R 

Mole fraction of oxidizer vapor at surface of oxidizer 
droplet . 

Array of values at each pressure must be entered for 
each temperature . 

Enter NPTP(l)xNPTT(l) number of values, 

6 per card with no embedded blank fields. 

Oxidizer heat of vaporization. Units: BTU/lbm 
Multiple arrays using same order as for cards 50, etc. 

Redlich-Kwong parameter "a" array used in the equation 
of state for oxidizer. Units: ft 4 -R^/lbm. 

Redlich-Kwong parameter, "b" array used in the equation 
of state for oxidizer. Units: fttylbm 


Cards 570-737: input oxidizer liquid heat capac- 

ity and oxidizer liquid enthalpy 


Number of pressures in liquid oxidizer table. 

Limit: 2 to 20 

Required to be zjro(O) 

Number of temperatures in oxidizer table. 

Limit: 2 to 20 
Required to be zero(O) 

Pressure array in ascending order for oxidizer table 
Units: psia 
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TABLF 

1 . 

INSTRUCTIONS FOR PROPF LEANT AND 
COMBUSTION CAS INPUT DATA (Cont.) 

: n. 

.**; .AT 

V/T. TABLE 
CCI -li 


DESCRIPTION 

590, 591, etc. 
I6E12.8) 

TTCPL(K.l) 
K=1 , NTCT ( 1 ) 


Tempernture in ascending order for oxidizer table. 
Units: °R 

600, 601, etc. 
(6F.12.8) 

(TCP(I ,K, 1) 
K=1 ,NTCP(1)) , 
I=-l ,NPCP(1) 

Specific heat at constant pressure of liquid oxidizer. 

Units: BTU/lbm-R. An array of values corresponding 
with temperature array must be entered to correspond 
with each pressure. Enter NPCP ( 1 ) xNTCP ( 1 ) number of 
values, 6 per card. Do not skip any fields. 

670, 671, etc. 
f6F12 . 8) 

(TH0L ( I , K , 1 ) 

K=1,NTCP(1)), 

T-1,NPCP(1) 

Enthalpy of liquid oxidizer. Units: BTU/lbm 

Multiple arrays using same input order as cards 600, etc. 



CarHs 940-1477: input tables of oxidizer and fuel 

vapor properties in which values of three dependent 
variables correspond with the same temperature array 
and at various pressure levels. 

940 

16112) 

NPV(l) 
NPV(2) 
NTV f 1 ) 
NTV (2) 


Number of pressures for oxidizer vapor tables. 
Limit: 2 to 20. 

Number of pressures for fuel vapor table. 

Limit: 0 to 20 

Number cf temperatures in oxidizer vapor table. 
Limit: 2 to 20 

Number of temperatures in fuel vapor table. 
Limit: 0 to 20 

950, 951, etc 
(6E12.8) 

TPV(K, 1) 
K=1 ,NPV (1) 


Pressure array in ascending order for oxidizer vapor 
table. Units: psia 

960, 961, etc 
(6E12.8) 

TTV(K, 1 ) 

K» 1 , NTV ( 1 ) 


Temperature array in asccgding order for oxidizer 
vapor tab 1 1 ?. Units: R 

970, 976, etc 
(6F13.8) 

(TCPV(K,I,1) 
1*1 »NTV(1)) , 
K*1 ,NPV(1) 

Specific heat at constant pressure for oxidizer vapor. 
Units: BTU/lbm. An array of values corresponding 

with temperature array must be entered to corres- 
pond with each pressure level. Do not ship any 
fields. 

1040,1041 

etc. 

(6E12.8) 

(TMUV (K, 1 , 1) 
I»1 ,NfV(l )) , 
K-l.NPV(l) 

Viscosity of oxidizer vapor. Uiits: lbm/ft-sec. 

Multiple arrays using same input order as for 
cards P?0, etc. 

1110,1111 

etc. 

(6EI2.8) 

(TH0V(K,I ,1) 
i-i.Nrv(i)), 

K-1,NPV(1) 

Enthalpy of oxidizer vapor. Units: BTU/lbm. 

Multiple arrays using same input order as for 
cards 970, etc. 




Omit cards 1250-1477 if NTV(2) - *. 


ri 


TABLE 1 . INSTRUCTIONS FOR PROPELLANT AND 

COMBUSTION GAS INPUT DATA (Cont.J 


VARIABLE 

CODE 



TPV(K,2) Pressure array in ascending order for fuel vapor table. 

K=1,NPV(2) Units: psia 


TTV (K, 2) 
K=1 ,NTV (2) 


Temperature array in ascending order for fuel vapor 


table. Units: R 


(TCPV (K, I , 2) 
1=1 ,NTV(2) ) , 
K=1 ,NPV (2) 

(TMUV (K, I , 2) 
1=1 ,NTV (2) ) , 
K=1 ,NPV (2) 

(T1I0V (K, 1 ,2) 
I=1,NTV(2)), 
K=1 , 1 ,NPV (2) 


NTDF(l) 

NTDF(2) 

TTDIF(I,1) 
1=1 ,NTDF(1) 


Specific heat at constant pressure for fuel vapor. 

Units: BTU/lbm. Multiple arrays using input order 

per card 970, etc. 

Viscosity of fuel vapor. Units: BTU/lbm. 

Multiple arrays using input order per card 970, etc . 


Enthalpy of fuel vapor. nits: BTU/lbm. 

Multiple arrays using input order per card 970, etc. 


Cards 1500-1554: input oxidizer tables of dif- 

fusion parameters. 

Number of temperatures in the oxidizer table. 

Limit: 2 to 20. 

Required to be zero(O) 

Temperature array in ascending order for oxidizer tabic 
Units: °R 


TDI FF ( I ,1,1), 
1=1, NTDF ( 1 ) 


TDIFF (I ,2,1) 
1=1, NTDF (1) 


Oxidizer binary diffusion parameter array (see page 86 ) 
for specie to stoichiometric combustion products. 
Units: ft 2 /sec 

Oxidizer binary diffusion parameter array for specie 
to fuel. Units: ft 2 /sec 


TDIFF (1,3,1) Oxidizer binary diffusion parameter array for specie 
1=1 ,NTDF (!) to oxidizer. Units: ft 2 /sec. 


(TPRF(1,K) Reference pressures used with the three corresponding | 
K*1 , 3) , oxidizer binary diffusion parameters. Units: psia 

(TTRF(1,K) Reference temperatures used with the three correspond- 

ing oxidizer bindary diffusion parameters. 'Jni*s:*R 





TAP IT. 1 


INSTRUCTIONS POP. FROITH-ANT AIT.' 

cm-wisi ion gas innn data (Cont.) 


\V:. AO. 

TO...V1' 


VAf. ' AivLP 
CCuii 


1610 

(6F12.8) 


1620 

(6E12.8) 


1630 

16F.12.fi) 


description 


rCRIT(l) 
TCP I . (1) 
PMWl.fl) 

F.MWV'(l) 


PCRIT (2) 
TCRIT (2) 
FMWL(2) 
F.MWV ( 2 ) 

ST0CMR 

FMWl’R 


Critical pressure of oxidizer. Units: psia. 

Critical temperature of oxidizer. Units: R 

Molecular weight of oxidizer as liquid. Units. Ibm/ 

lb-mole . . 

Molecular weight of oxidizer as vapor. Units: Ibm/ 

1 b mo 1 e 

Critical pressure of fuel. Units: psia 

Critical temperature of fuel. Units: R 

Molecular weight of fuel as liquid. Unit^: ibm/lb-mole 

Molecular weight of fuel as vapor. Units: Ibm/lb-moK 

Stoichiometric mixture ratio 
Molecular weight of products at. ST0CMR . 

Units: Ibm/ lb -mole 


1640 

(6112) 

1650, 1651 , 
etc. (6H12.fi) 

1660,1661 
etc. (6L12.8) 

1670,1671 
etc. (6F12.8) 


1680,1681 
etc. (6F.12.8) 


1690,1691 
etc. 

(6F12.8) 

1700,1701 
etc. (6E12.8) 

I 
I 


NMRCf.F 

NTCGF 

TMRCGF(l) 

1 = 1 .NMRCC.F 

TTCCP(I) 

1=1 .NTCGF 

tmwcgf (1 , j) , 

.1=1 .NTCGF 


TMUCGF (1 ,.l) 
J=1 .NTCGF 


TCPCGF (1 , J) 
J=l, NTCGF 


TMWCGF(2,J), 
J-l, NTCGF 

I 

I 


Cards 1640-2254: input table of combustion gas 

properties with mixture ratio and temperature as 
tne independent variables. This table is used to 
determine droplet film properties at the mean 
temperature between the droplet and free stream gas 

Number of mixture ratio levels Limit: 2 to 20 

Number of temperatures at each mixture ratio level. 
Limit: 2 t > 20 

Mixture ratio array in ascending order. 


Temperature array in ascending order. Units: 


Molecular weight array for combustion gas at the first 
mixture ratio level and corresponding with TTCGF array. 
Units: ibm/lb-mole 

Viscosity array for combustion gas at the first mixture 
ratio level and corresponding with TTCGF array. 

Units: lbm/ft-sec. 

Specific heat array for combustion gas at the first 
mixture ratio level and corresponding in order with 
TTCGF. Units: BTU/lbm-R 

Repeat arrays of TMWCGF, TMWCGF and TMUCGF for each 
mixture 


TABLE 1 . INSTRUCTIONS FOR PROPELLANT AND 

COM'USTION GAS INPUT DATA (Concluded) 


CARD NO . 
r, FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

1 

1 

etc. 

i 

etc . 

Unit conversion options: 



If TMUCGF (1,1) is negative then the TMUCGF array is 



divided by 3600. 



If TMUCGF (1,2) is negative, then the TMUCGF array is 



multiplied by 32.16. 



r 

Cards 2260-2427: input tables of oxidizer liquid 



surface tension and viscosity as functions of tempera- 



ture and pressure. 

2260 

NPST(l) 

Number of pressures in oxidizer table. Limit: 2 to 20. 

(6112) 

NPST (2) 

Required to be zero (O) 


NTST(l) 

Number of temperatures in oxidizer table. Limit: 2 to 20 


NTST(2) 

Required to be zero (0) 

2270,2271 

TPST(K, 1) 

Pressure array in ascending order for liquid oxidizer. 

etc. (6E12.8) 

K=1 ,NPST(1) 

Units: lbf/in. 2 . 

2280,2281 

TTST (K, 1) 

Temperature array in ascending order for liquid. 

etc. (6E12.8) 

K=1,NTST(1) 

Units: R. 

2290,2291 

(TST (I , K, 1) 

Oxidizer liquid surface tension. Units: lbf/ft. 

etc. (6E12.8) 

K=1 ,NTST(1) ) , 

An array of values corresponding with temperature array 


1=1 ,NPST(i ) 

must be entered to correspond with each pressure level. 



Do not skip fields. 

2360,2361 

(TVISL (I , K, 1) 

Oxidizer liquid viscosity. Units: lbm/ft-sec 

etc. (6E12.8) 

K= 1 , NTST ( 1 ) 

J > 

Multiple arrays using same input order as cards 2290, etc. 


1=1 ,NPST(1) 




Cards 2600-2687: input tables of fuel compressi- 



bility factor as a function n f temperature and pressure. 

2490 

NPZ (] ) 

Required to be zero (0) 

(6112) 

NPZ (2) 

Number of pressures in fuel table. Limit: 2 to 20. 


NTZ(l) 

Required to be zero (0) 


NTZ(2) 

Number of temperatures in fuel table. Limit: 2 to 20. 

2600,2601 

TPZ (K, 2) , 

Pressur. array in ascending order for fuel. 

etc. (6E12.8) 

K=1,NPZ(2) 

Units: psia 

2610,2611 

TTZ (K, 2) 

Temperature values in ascending order for fuel. 

etc. (6E12. 8) 

K=1 ,NTZ (2) 

Units: R 

2620,2621 

(TZ (I , K, 2) 

Compressibility factor of fuel. An array of values cor- 

etc. (6E12.8) 

K=1 ,NTZ (2) ) 

, 

ponding with temperature array must be entered to cor- 


1=1 , NPZ (2) 

pond with each pressure value. Do not skip fields. 
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TABLE 



1 


INSTRUCTIONS FOR CONTROL INPUT DATA 



CARD NO. 


VARIABLE 

CODE 



10 


30 

(6E12.8) 


40 

(6E12.8) 


SO 

(6E12.8) 


ICUPC 


NCHAMC 


NC0N4C 


EMRGJC 

STGJC 

EMWGJC 

GAMGJC 

XLMC 


DELTXC 


BSPRC 


CSPRC 


XMINDE 


CLNTC 

C0NRAC 

CCANGC 


Control indicator: value <_ 0 to bypass DER option; 

value > 0 for number of injector flow zones analyzed 

Control indicator: value of "0" for one cup or a 

chamber calculation, value of "1" for both cup and 
chamber calculations. If IDER > 0, program sets 
ICUPC = i. 

NOTE: If IDER = 0 and ICUPC = 0, then NCHAMC, M2C and 

NC0N4C are ignored. 

Control indicator for type of chamber geometry input: 
value of "0" for conventional geometry (card 50), 
or an integer for the size of a cross-sectional 
area array (card 60, 61, etc.) 

Print control indicator: solution printed at calculation 

step intervals of M2C. 

Print control indicator to force printout of each step 
for first NC0N4C chamber ca lculations . ___ 

Include cards 30-60 only if either IDER or ICUPC > 0. 

This card group is used to define chamber parameters 

when the computer run includes the analyses of both 

injector element cup(s) and chamber. 

Total chamber "rigimesh" (or gas mantle) flowrate at 
injector face. Units: lbm/sec. 

Weight mixture ratio (oxidizer/ fuel) of WGJC flow. 

Stagnation temperature of WGJC. Units: R. 

Molecular weight of WGJC. Units: lbm/lb-mole. 

Specific heat ratio, Y, of WGJC. 

Length of mixing region. Rigimesh flow is mixed into 
the combustion region linearly over this region. 
Units: in. 

Axial step size for chamber calculations. Units: in. 

Recommended value = 0.05 in. 

Droplet formation size parameter in chamber.* 

Recommended value = 120.0 

Liquid jet stripping rate parameter in chamber.* 
Recommended value = 0.08 

Minimum axial distance for DER punch card output (not 
required if IDER <_ 0). Units: in. 

Include card 50 only if NCHAMC - 0. 

Cross-sectional area of chamber at the injector end. 
Units: sq. in. 

Chamber length from injector face to the throat plane. 
Units: in. 

Chamber contraction ratio (area of c'namber/area throat) . 

Nozzle angle of convergence. Units: degrees 


‘See text page 


for description. 


i 





. INSTRUCTIONS FOR CONTROL TNPUT DATA fConcT U 


50 

fCont . ) 


60, 61, etc. 

(6E12.8) XCUAMC (1) 





TABLE 4 . INSTRUCTIONS FOR CASE INPUT DATA 


V.. \ . , 1 ) ,\ o . 



no 

(18A4) 

111 

(18A4) 

120 

(6112) 


130 

(6112) 


VAR I AD 
CODE 


AMAT(I) 

1=1,18 

AMAT(I) 

1=19,36 

NDSCI 

NELEM 

NCUAM 


ICUP 


ICPE 


I READ 


M2' 

NC0N4 


DESCRIPTION 


At least one group is always required. 

If ID) ,R > 0, then IDER number of groups are included, 
each one of which define injector element param- 
eters for a specific flow zone. 

If IDER <_ 0 and ICUPC = 0, then only one of these groups 
are included which defines parameters either for a 
cup analysis or for a chamber analysis. 

If IDER 0 and ICUPC ~ .1, then only one of these groups 
are included which defines parameters for a cup 

analysi s . 


Case comment card one. 


Case comment card two 


Number of spray drop sizes at start plane. 

Number of injector elements in case. 

Control indicator for type of injector cup (or chamber) 
geometry input: value of "6" for conventional 

geometry (card 140) or an integer for the size of 
a cross -sectional area array (Card ISO, 151, etc.) 
Control indicator on type of case analysis: 

"1" for an injector cup analysis only 
"2" for both an injector cup and a chamber analysis 
"3" for a chamber analysis only 
Control indicator on gas expansion: 

"1" for constant pressure expansion limitation 
"0" for lull expansion to fill cross-sectional area 
Control indicator which specifies input to be read 
for next case: 

"0" to start with control card 10 
"1" to start with case card 110 
”2" to start with combustion gas (CGTBIN) table 
input 

"3" to start at beginning of input (includes 
TABIN and CGTBIN tables) . 

Print control indicator solution printed at calculation 
step intervals of M2. 

Print control indicator to force printout of each step 
for first NC0N4 axial calculations. 



TABLE 4 . INSTRUCTIONS TOR CAST INPUT DATA fCont.l 


('A!. U .'0. 

,, format 

VARIABLE 

CODE 

DESCRIPTION 

1^0 fcont.*) 

TFYrGI, 

Control indicator for expansion around liquid post ffor 
cup calculation only)* 

"1" constant gas expansion 

"2" liquid and gas expansion 

"3" liquid expansion and gas contraction 

Recommended: IEXPGL = 3 , 


IAT0 

Atomization control indicator. Fnter value of "1". 

140 


Include card 140 only if NCI JAM = 0. 

(6F.12.8) 

ACS I 

Cross-sectional area of injector cup or chamber at 
upstream end. Units: sq. in. 



Injector cup or chamber length. Units: in. 


1 

Area ratio of injector cup or chamber: ACSI over cup 

exit or nozzle throat area. 


CCANG 

Angle of convergence: for cup, a negative value 

specifies angle of divergence; for chamber, value 
is nozzle angle of convergence. Units: degrees. 


RCBC 

Wall radius of curvature leading into convergent section. 
Units; in. For cup, set RCBC = 0. 


RCT 

Wall radius of curvature entering nozzle throat. 
Units: in. For cup, set RCT = 0. 

ISO, 151, etc. 


Include th se cards only if NCUAM >0. 

(6F.12.8) 

XCHAM(l) 

First value in array of axial distances from the begin- 
ning of either injector cup for cup analysis or 
iniector face for chamber analysis. Units: in. 


ACHAM(l) 

First value in array of cross-sectional areas corres- 
ponding with position XCHAM(l). Units: sq. in. 


XCIIAM(2) 

Fnter NCUAM pairs of values, 3 per card, with XCIIAM in 


ACUAM(2) 

1 

1 

etc. 

ascending order. 

160 

WCGI 

Flowrate per element of gas stream surrounding liquid 

(6E12.8) 


jet at start position of case. Units: lbm/sec. 


FMRCGI 

Weight mixture ratio (0/F) of WCGI. 


ACGI 

Initial cross-sectional flow area of WCGI. Units: sq.in. 


EMRII 

Weight mixture ratio of gas in manifold. EMRII = FMRCGI 
when gas is fully reacted at start position of case. 


STT 

Stagnation temperature of WCGI at a reference mixture 
ratio AMRT. Units: °R 


AMRT 

Reference weight mixture ratio for WCGI temperature of 
STT. 


♦See text page 22 for descr r , ion. 
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1 



TABLE 4 . INSTRUCTIONS FOR CASE INPUT DATA (Cont.) 


CARD NO. 
£• FORMAT 

VARIABLE 

CODE 

DESCRIPTION 

170 

WLJI 

Flowrate per element of oxidizer liquid jet at start 

(6E12.8) 


position. Units: lbm/seg. 

Temperature of WUI. Units: R 


TLI 


VLJI 

Velocity of WLJI. Units: ft/sec. (If VLJI < 0, area 

of WLJI. Units: sq in.) 


D0DMAX 

Maximum dropsize permitted in atomization of liquid jet. 
Units: microns. 


BSPR 

Droplet formation size parameter.* 

Recommended: BSPR = 120.0 (Chamber) ,3.0553 (Cup) 


CSPR 

Liquid jet stripping rate parameter.* 

Recommended: CSPR = 0.08. (Chamber), 0.037854 (Cup) 

x80 

WGJI 

Flowrate per element of gas stream surrounding WCGI at 

(6E12.8) 


start position. Units: lbm/sec. 


EMRGJ1 

Weight mixture ratio (0/F) of WGJI. 


STGJ 

Stagnation temperature of WGJI. Units: R 


EMWGJ1 

Molecular weight of WGJI. Units: lbm/lb-mole. 


GAMGJ1 

Specific heat ratio, y, of WGJI. 


XLM 

Length of region for WGJI to be mixed with WCGI. 
Units: in. 

190 

PCI 

Injector end static pressure. Units: lbf/sq.in. 

(6E12.8) 

CUPDP 

Estimated static pressure drop in injector cup. For 
cup analysis only. Units: nsi 


CUPDPL 

Tolerance on matching cup exit pressure with PCI. 
Units: psi. 


STX2 

Start plane position; either distance from liquid 
injection post for cup or distance from 
| injector face for chamber. Units: in. 


DELTX2 

Axial step size for case computations. Units: in. 


FCHA 

Fraction of chamber cross-sectional area taken by this 
flow zone case. 

191 

RFLAME 

Radial location of the pseudo flame front. Recommended: 

(6E12.8) 


RFLAME = radius of the fuel sleeve. Units: in. 


XFLAME 

Axial location of the start of the pseudo flame front. 
Recommended: XFLAME = 0.0 (injector face). 

Units: in. 


VFLAME 

Turbulent flame speed. Units: ft/sec 

200,201 


Include these cards only if NDSCI > 0 

etc. (4E12.S; 

V0DI (1) 

Droplet velocity of spray group 1. 
Units: ft/sec 


T0DI (1) 

Droplet temperature of spray group 1. 
Units: R 


D0DI(1) 

Droplet diameter of spray group 1. 
Units: microns 


WSPRI(l) 

Spray group 1 flowrate. Units: lbm/sec 

Enter NDSCI number of spray groups. 


*See text page for description. 




INSTRUCTIONS FOR CASE INPUT DATA (Concluded) 


DESCRIPTION 


Include cards 300-331 only if IDER > 0 

Number of mixing zones per element (maximum of 40) . 

Maximum number of oxidizer droplet spray groups for 
DER punched output (maximum of 11). 

Fraction of total case fuel flowrate in the first mixing 
zone. 

Fraction of total case oxidizer flowrate in the first 
mixing zone. 

Enter NMIXZ pairs of values, 3 per card. 

Fraction of total spray flowrate in the first DER spray 
group. 

Enter NG0 values, 6 per card. 





6E12.8 


Decimal variables with variable names beginning 
with letters other than I through M 
(Use decimal point or account for implied decimal 
location, one value every 12 spaces, 6 consecutive 
values per card.) 

The "VARIABLE CODE" column gives the FORTRAN code names of input variables 
as they appear in the program listing. A single value is to be entered 
for each coded variable unless it is subscripted. Array sizes for sub- 
scripted integer and decimal variables are also indicated within paren- 
thesis in this column, following the variable name. For most of the data, 
all of the values of one variable are read before proceeding to the next 
variable. Note that some arrays with multiple subscripts are "packed", 
i.e., values for each subscript level start immediately in the next field, 
not skipping fields to start on a new card. 


Variable names and/or descriptions of variables are given together with 
appropriate dimensions and limits, in the "DESCRIPTION" column. Generally, 

the program is written in units of lb-in-sec-°R, but there are some 
exceptions. 


PROPELLANT AND COMBUSTION GAS INPUT DATA 

The first block of data required as input to the CICM computer program 
comprises the propellant find combustion gas properties (Table 1). Print 
out of tnis block during execution is controlled by the variable IPTAB: 
"0" to suppress printout, "1" to print the data block. 
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Liquid, Vapor and State Properties 
of Propellants (Card No. 20 et seq,._) 

F.xtensive tables of propellant properties are provided as input to t! ‘ 
droplet diffusion model. Tbe first of these tables gives values for 
the vapor mass fraction, X y , at the droplet surface (equivalent to a 
reduced partial pressure), the heat of vaporization, AH y , and parameters 
a and b of the Redlich-Kwong equation of state. Tables of X y and AH y 
as functions of both total pressure and temperature are provided, while 
a and b are provided as functions of temperature only. As noted in 
Table 1 , only oxidizer properties are required. 

Values in these tables should correspond to temperatures ranging from 
injection temperature to the critical temperature only. Pressure ranges 
should cover the pressure variation occurring in the subsonic flow por- 
tion of a combustor under analysis. It is probably preferable to input 
data for much wider variation so that the same tables can be used for 
other engines using the same oxidizer. This approach was taken in struc- 
turing the liquid oxygen tables supplied with the example case, Appendix B. 

Values of X y and AH y should include real gas effects, i.e., dependence 
upon total pressure level. For vapor-liquid equilibrium, the free energies 
of the vapor and the liquid are equal. This fundamental relationship for 
vapor- liqui I equilibrium is conveniently expressed in terms of fugacities; 
for each component i the fugacity of the vapor, f.\ is equal to that of 
the liquid, f^, (Ref. 7 ). Because the liquid senses the total pressure 


R3 

_ :_l_j 


while the vapor senses only its partial pressure, the equilibrium 
relationship may be written as 


*1 ( p v.’ ■ f i l (P Totai> 

Hence, at constant temperature, as the total pressure increases the 
partial pressure of the equilibrium vapor also increases. 


In the calculation of vapor-liquid equilibrium, the vapor must be con- 
sidered a non-ideal gas. Of the four two-constant equations of state 
which have been widely used, the Redlich-Kwong equation is accurate 
throughout the pressure and temperature range and is the most accurate 
at high pressures. The Redlich-Kwong equation is: 

p P-T a 

" (v-b) T 0.S y(v+b) 

The parameters a and b are determined from mixing rules (Ref. 7 ). To 
match data over wide ranges, a and b may be expressed as functions of 
temperature. 

Data for these tables may be obtained by solving simultaneously four 
equations given in Ref. 7 > which are expressions for the liqu'd and 
vapor fugaoities and liquid and vapor spates. Note that, at super- 
critical pressures, AH y •+• 0.0 at temperatures well below the critical 
temperature . 


84 


For a non-ideal gar the species vapor enthalpy is a function of its 
partial pressure in the gas (Ref. 10), and is thus depen lent on the 
total pressure. Hence, the heat of vaporization 

v vapor liquid 

is a function of total pressure as well as cf liquid temperature. 

Liquid Specific Heat and 
Enthalpy (Ca rd No. 570 et seq.) 

The next olock of propellant property data provides liquid specific 
heat and liuuid enthalpy as functions of pressure and liquid tempera- 
ture. Again, only oxidizer properties are required. Note that, 
although these are denoted as "liquid" properties, the tables should 
provide data to temperatures as high as the combustion gas temperature; 
for temperatures higher than the saturation temperature corresponding 
to the tabulated pressure, the pure vapor properties are used. 

Vapor Specific Heat, Viscosity, an d 
Enthalpy (Card No. 940 et seq.) 

The next data to be input are tables of vapor specific heat at constant 
pressure, vapor viscosity, and vapoi enthalpy as functions of pressure 
and temperature. These may be derived from tabulations of experimental 
data or from standard correlation methods, c.g., such as those given in 


R5 


Ref. 10 . 


T 


Binary Diffusion Coefficient 
Parameter (Card No. 1500 eg seq.) 

Oxidizer binary molecular diffusion coefficients are calculated in the 
program from the data input in the TDIFF (I, K, 1) tables. This param- 
eter is assumed to be a function of temperature; tabulated values cor- 
respond to temperatures in the array TTDIF(I , 1). The subscript I 
denotes the various temperature levels. The subscript V indicates the 
gaseous specie into which the oxidizer is diffusing into, as noted in 
the description in Table l . 


The TDIFF parameter has the following meaning: An equation for binary 
diffusion coefficients, based on use of the Lennard-Jones potential in 
a kinetic theory model, is given in Ref. 10 as: 


D 


2 

12 


0.001858 



♦ M 2 )/Mj m 2 



Multiplying and dividing this equation by a reference temperature and 
reference pressure gives: 


0.001858 T 


3/2 


°12 


ref 


P ref °12 


[(Mj ♦ M 2 )/M 1 M 2 1 



P re*\ 


The product 


0.001858 T r#f 3/2 [(Mj ♦ M 2 )/Mj H ? 


ref *12 


is tabulated as the TDIFF parameter. 
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This is assumed to vary with temperature, but not with pressure. 


Note that TDTFF (I, 1, 1 ) are for the oxidizer species diffusing 
into combustion products at stoichiometric mixture ratio. For 
lower or higher mixture ratio combustion gases, the multicomponent 
diffusion coefficient is approximated by the program for the oxi- 
dizer species diffusing into a mixture of stoichiometric products 
and excess fuel or oxidizer vapor, respectively. 

Propellant Critical Properties and 
Molecular Weight (Cards No. 1610, 1620, 1630) 

The critical temperature, critical pressure, stoichiometric mixture 
ratio, and molecular weight of the stoichiometric products are input 
in this data block. The vapor molecular weight will differ from that 
for the liquid only if there is vapor phase decomposition. If this 
occurs, it is recommended that the heat cf dissociation be included 
in the tabulated values of heat of vaporization. 

Combustion Gas Properties at 

Film Conditions (Card No. 1640 eq seq.) 

Combustion gas film properties required in subroutine FGTR0P for cal- 
culating film gas properties are molecular weight, viscosity, and 
specific heat. These are tabulated as functions of mixture ratio and 
gas temperature. For the oxygen/hydrogen data deck supplied with the 
sample case, these data were obtained from the Rocketdyne free energy 
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equilibrium performance program by specifying different values of 
mixture ratio and product temperature (rather than mixture ratio and 
initial enthalpy). 

Oxidizer Liquid Surface Tension 
and Viscosity (Card No. 2260 et seq.) 

The next propellant properties to be input are tables of liquid surface 
tension and liquid viscosity as functions of pressure and temperature. 

The tables should include temperatures ranging from injection tempera- 
ture to the oxidizer critical temperature. 

Fuel Compressibility Factor 
(Card No. 2490 et seq.) 

Tables of fuel compressibility factor are input as a function of pres- 
sure and temperature. The tables should include temperatures ranging 
from fuel injection temperature to the combustion gas temperature. 

STAGNATION EQUILIBRIUM COMBUSTION GAS INPUT DATA 

The second block of data required as input to the CICM computer program 
comprises the stagnation equilibrium combustion gas (Table 2). Com- 
bustion gas properties, tabulated as functions of gas mixture ratio, are 
obtained from prior peripheral computation using a thermodynamic equilib- 
rium performance program. Rocketdyne's free energy performance program 
was used to generate the table supplied in the reference case, but any 
comparable program would be sufficient. The combustion temperature. 
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molecular weight, specific heat, and viscosity entered in this table 
are properties for equilibrium combustion products at stagnation con- 
ditions corresponding to the mean expected chamber pressure. The 
properties are assured to be functions only of mixture ratio and not 
pressure. 


CONTROL INPUT DATA 


The third block of data required as input to the CICM computer program 
comprises control data and also chamber conditions (Table 3 ). 


Indicator Card (Card No. 10) 

The first control input data card contains indicators for controlling: 
(1) the DER option (IDF.R), (2) coupled cup-chamber calculations (ICUPC), 
(3) the type of chamber geometry input (NCHAMC) , and (4) the chamber 
solution printout intervals (M2C and NC0N4C) . p or execution of the 
program using the DER option, IDER specifies the number of injector flow 
zones (or number of different element types) to be used in the analysis 
of the engine. If IDER and ICUPC are both less than or equal to zero, 
this card is the only control card required as input. 


"Rigimesh" (or Gas Mantle) Conditions 


(Card No. 30) 

The next card of control data specifies the "Rigimesh" flow conditions in 
the chamber. Even if the "Rigimesh" flowrate is zero, it is recommended 
that values for the "Rigimesh" stagnation temperature (STGJC), molecular 
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weight (EMWGJC), and specific heat ratio (GAMGJC) be entered in order 
to avoid possible execution errors. At present, the CICM program 
mixes the "Rigimesh" flow into the combustion region linearly with 
position from the injector face to the axial location specified by 
XLMC . 


Parameter Card (Card No. 40) 

The next control data card specifies the axial step size for the cham- 
ber calculations (DELTXC) > the chamber droplet formation size and liquid 
jet stripping rate parameters (BSPRC and CSPRC) , and the minimum axial 
distance for DER punched card output (XMINDE) . 

During execution of the program using the DER option, all IDER flow 
zones are executed to an axial location specified by the length of the 
longest liquid jet or to the axial distance specified by XMINDE, depend- 
ing on which is larger, before DER punched cards are generated. 

Chamber Geometry (Card No. 50, 60, etc.) 

The last set of control cards specify the chamber area as a function of 
axial distance. Two different methods of input are possible, depending 
on the value of NCHAMC. The first method, NCHAMC<0, requires the 
cross-sectional area of the injector face (ACSC) , chamber length (CLNTC) , 
chamber contraction ratio (C0NRAC) , nozzle angle of convergence (CCANGC) , 
and the radii of curvature of the beginning of convergence (RCBCC) and 
at the throat (RCTC) to describe the combustor area as a function of 
length. 
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With the second method, NCHAMC>0, the geometry of the combustor is 
specified through the array ACHAMC. At selected axial positions 
(XCHAMC), the chamber e~ea is given by AC11AMC. For axial locations 
between the selected values of XCHAMC, the program linearly inter- 
polates for the combustor area. 

CASE INPUT DATA 


The final block of data required as input to the CICM computer program 
comprises the case input data (Table 4 ). Tf the DF.R option of the pro- 
gram is being utilized, the program requires TPFR number of case-input 
data blocks. 

Comment Cards (Card No. 110 and ill) 

Two alphanumeric (A -formatted') comment cards are provided to permit the 
user to document the case with such information as injector name, draw- 
ing number, element description, propellant combination, nominal chamber 
pressure and mixture ratio, date of the computer run, etc. 

Indicator Cards (Card No. 120 and 130) 

The next two case data cards contain variables (indicators) for: 

(1) specifying the number of spray drop sizes at the start plane (NDSCT) , 

(2) specifying the number of injector elements in the case (NELEM) , 

(3) controlling the type of injector cup (or chamber) geometry input 
(NCHAM) , (4) controlling the type of case analysis (ICUP) , (4) controlling 
the gas expansion for the first incremental step (ICPE), (5) controlling 
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the input to be read for next case (IREAD), (6) controlling the case 
solution printout intervals (M2 and NC0N4) , (7) specifying the method 
of describing the gas expansion around the liquid post (IEXPGL), and 
(8) controlling the atomization process (IAT0) . A constant pressure 
expansion option (ICPE = 1) has been included in the CICM computer 
program to allow the combustion gas, in the absence of "Rigimesh" flow, 
to expand at constant pressure in the chamber. For cup calculations, 
this indicator should be set equal to zero. The different options for 
expansion around the liquid oxidizer post (IEXPGL) are discussed on 
page 22. This indicator is required only for cup calculations. 

Case Geometry (Card No. 140. 150, etc.) 

The next set of case cards describe the flow area as a function of 
axial distance. This set of cards is very similar to the chamber 
geometry cards described in the control input data block. For most co- 
axial engines, the case geometry will describe the cup geometry. Two 
different methods of input are possible depending on the value of NCHAM. 
The first method, NCHAM <0, uses the cross-sectional area of the 
injector cup or chamber at the upstream end (ACSI) , injector cup or cham- 
ber length (CLNT) , ratio of inlet area to exit area of the injector 
cup or chamber (C0NRAT) , angle of convergence (negative value specifies 
angle of divergence), and the radii of curvature at the beginning of 
convergence (RCBC) and at the injector cm or chamber exit (RCT) . 
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In the second method, NCHAM>0, the geometry of the injector cup or cham- 
ber is described through the array ACHAM. At selected axial positions 
(XCHAM), the injector cup or chamber cross-sectional area is ,»iven by 
ACIIAM. For axial locations between the selected values of XCHAM, the 
program linearly interpolates for the injector cup or chamber area. For 
cup calculations, XCHAM is the distance from the upstream end of the 
injector cup. For chamber calculations, XCHAM is the distance from the 
injector face. 

Combustion Gas Conditions (Card No. 160) 

The next card for the case data specifies the combustion gas, or fuel , 
flow conditions at the computational start plane. The combustion gas, 
or fuel, flowrate at the start position (WCGI) , the weight mixture ratio 
at the scart plane of the gas (FMRCGI) , the initial cross-sectional flow 
area of the gas (ACGI), the weight mixture ratio of the gas in the mani- 
fold (FMRII), and a reference stagnation temperature (STT) and mixture 
ratio (AMRT) . For cup calculations, the mixture ratio at the start 
plane (F.MRCGT) will be the same as the mixture ratio in the manifold 
(FMRII) . 

For chamber calculations, generally the mixture ratio at the start plane 
will not be equal to the manifold mixture ratio. The reference tempera- 
ture (SIT) and reference mixture ratio (AMRT) are used to update the 
equilibrium stagnation gas tables (Table 2 ) to account for differences 
in propellant energies. Nominally, these reference values are equal to 
the combustion gas, or fuel, manifold stagnation temperature and mixture 
ratio. 


93 


T 


T 


Liquid Jet Conditions (C ard No. 170) 

Thc next card for the case data specifies the liquid Jet flow conditions 
the start plane and droplet parameters. The liquid Jot flowrate 

temperature (TLI) , velocity (VLJI), maximum dropsize permitted 
in atomization of the liquid jet (D0DMAX) , and parameters describing the 
droplet formation size (BSPR) . and liquid jet stripping rate (CSPR) are 
included. If a negative value is input for the liquid jet velocity. 

(VLJI), the program will interpret this value to be the liquid jet cross- 
sectional area CALJI . -VLJI). If the local droplet diameter produced 
by the stripping process is larger than D0DMAX, stripping of the liquid jet 
will cease until the local droplet diameter is smaller than D0DMAX. 


"Rjpimesh" (or Gas Mangle) Condit ions, 

(Card No. ISO) 

Th0 next case data card specifies the "Rigimesh" flow conditions in the 
chamber This card is very similar to the ••Rigimesh" condition card 
described in the control input data bloch. Although for cup calculations, 
the "Rigimesh" flowrate (WGJI) must be set equal to zero, 

mended that arbitrary values for the stagnation temperature (STGJ). molec- 
lar weight (EMKGJ1) , and specific heat ratio (GAMGJ1) be entered to 
,void possible execution errors. At present, the CICM program mixes the 
"Rigimesh" flow into the combustion region linearly from the injector 
face to the axial location specified by XLM. 
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Pressure and Distance Card (Card No. 190) 


The next case data card specifies the injector face static pressure 

(pen, estimated rup static pressure drop fCIIPOP) . the tolerance oa 

matching the cup exit pressure with the injector face pressures 

(CUPDPL) , start plane position (STX2) , axial step size (0FI.TX2) , and 

the fraction of the chamber cross-sectional area represented by this 

flow zone case (FCHA) . For cup calculations, it is recommended that 

the start plane location be chosen as equal to one liquid post thick - 

na „ r< . T v 7 _ t ) and the axial step size be set equal to 

ness (STX2 - t liquid pos t ; 

0.00S inch. For chamber calculations, it is recommended that the 
stnrt plane he chosen a, the injector face fSTX? - 0.0) and the axial 
step size be set equal to 0.05 inch. 

F l am e Propagation Conditions (Car d No. 191) 

The next case data card specifies the radial location of the pseudo- 
flane front (RFLAMF.) , the axial location of the start of the flamefront 
relative to the injector face (XFLAME) , and the turbulent flame speed 
(VFLAME). It is recommended that the radial location of the .Jamefront 
be set equal to the fuel sleeve radius and the axial location of the 
start of the flamefront be set equal to zero (injector face). If mF.R 
and NDSCI are both less than or equal to zero, this card is the last 

case data card. 
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Droplet Spray Grot ip Dcscripti o n 
(Card No, 200, etc.) 

Ihe next set of case data cards specify the droplet spray groups present 
at the start plane. This set of cards is needed only if NDSCI>0. The 
droplet velocity (V0DI), temperature (T0DI) , diameter (D0DI) , and spray 
group flowrate (WSPRI) are included. Each droplet group is entered on 
a separate card, i.e., there will be NDSCI cards in the set. 

PER Parameters (Card No. 300, etc.) 

The last set of cards for the case input data specify the parameters used 
m interfacing the CICM program with the STC section of DER. These cards 
are included only if IDER>0. Included is the number of mixing zones 
per element (NMIXZ) , the number of oxidizer droplet spray groups for 
PER punched card output (NG0) , fractions of the total fuel and oxidizer 
flowrate for the case in each mixing zone (FFMIX and F0MIX) , and the 
fraction of total spray flowrate for the case in each DER spray group. 
Cold flow data are required to dofine the number of mixing zones and the 
fractions of fuel and oxidizer in each zone. This is the last set of 
cards in the case input data block. 


PROGRAM OUTPUT 



The output of the Cl CM computer program is provided as the usu ? 1 t-r-u’ar 
printout. A sample case is included in Appendix C. Input da ' prin ed 
as they are read which permits both a full documentation of thfc ror;:u t o' 
run conditions for later analysis and a convenient method to check the 
input for errors if unusual results are calculated. The input sections 
should be examined for each case run to be sure that the intended input 
data were actually used. 

During CICM analysis, data are written out as they are generated. At 
selected axial incremental positions, complete gas and propellant spray 
group data are printed. Additionally, the percentages of propellants 
atomized, vaporised and reacted are listed. At the top of each axial 
station printout, two comment cards (from the case data input block) are 
listed and an identification line is written to inform the user whether 
the calculation was a cup or chamber case. 

Upon completion of the case, the program writes out an identification line 
informing the user that the case calculation is finished. For calcula- 
tions, the program checks to see if the cup exit pressure is equal to the 
chamber pressure. If the two pressures are different, the program prints 
out this fact and reruns the cup case with a new estimated cup delta 
pressure. 
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If the DER option has been specified, after all zone cases have been 
executed the program checks to see if each chamber case was continued 
to the axial position required for DER punched card output. If any of 
the zone chamber cases was terminated before reaching the axial position 
required for DER punched card output, the program automatically recal- 
culates these zone chamber cases. Upon completion of all zone calcula- 
tions, the program lists the DER punched card output. 

Upon analyzing all of the input data, the program writes an identification 
line informing the user that the program is terminating in the normal 
fashion. 



f'RROR ANAl.YSJS 


The most common cause of errors during execution of the CICM computer 
program are mistakes in the input information. The program contains 
certain special printouts if input limits are exceeded (subroutine 
TABIN) and if interpolation beyond reasonable limits of tabulated 
tables is attempted (subroutine I-flfTAC) . The usual reason for these 
error messages is bad input data. 

If the model calculations are allowed to proceed to the nozzle throat, 
the program may terminate calculations before the throat plane is 
reached if the calculated combustion gas velocity exceeds the local 
sonic velocity. This early termination will not effect sequenced model 
calculations and should not be encountered during execution with tie HER 
option. The early termination can be corrected by adjusting the injector 
static pressure . 

In oxecuting the program with the DER option being used, the user should 
verify that the number of DER droplet spray groups (variable NG0) is th" 
same for all DF.R zones. If the number of spray groups arc different, 
the resulting punched card output will be inconsistent. 
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APPENDIX A 
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APPENDIX C 



CONTENTS OF 
SAMPLE CASE OUTPUT 


V** " .-fiw 


Approx. No. 
of Pages 

Title Page Identifying Computer Model 1 
Propellant and Combustion Gas Input Data 42 
Control Innut Data 1 
Cup Calculation of Element Type #1 


1st Axial Pass 

Case Input Data (cup) 2 

Solution at Selected Axial Stations 18 

2nd Axial Pass with Estimated Cr.p AP Adjusted 

Modified Case Input Data 2 

Solution at Selected Axial Stations 18 

3rd Axial Pass with Estimated Cup AP Adjusted 20 

Chamber Calculation for Single Element Type #1 

Case Input Data (chamber) 2 

Solution at Selected Axial Stations 18 

Cup Calculation of Element Type #2 

1st Axial Pass 10 

2nd Axial Pass 10 

Chamber Calculation for Single Element Type #2 20 

Renin of Chamber Calculation for Single Element Typo #1 to 

Extend Distance to Match That Required for Element Type #2. 

Type #2 element Liquid Jet Was Greater Than Distance 

Specified in Input (XMINDE) 22 

Gas and Spray Data Which Were Calculated 8 Punched for 

DER/STC Program 2 


C-ii 


SUMMARY OF SAMPLE CASE SPECIFICATION 


Injection Elements 
Number of Elements 

Distance From Convergence of Streams 
to Injector Face 

Cup Cross-Sectional Area 
Depth of Cup Flare 

Liquid Jet Injection Area 

Annular Injection Area of Gas Stream 

Liquid Jet Flowrate 

Gas Stream Flowrate 

Gas Stream Mixture Ratio 

Rigimesh Gas Flowrate 

Chamber 

Area at Injector Face 
Area at Throat 
Chamber Length 
Chamber Shape 


Units 

Type #1 

Type #2 

- 

30 

36 

in. 

0.10 

0.10 

. 2 
in. 

0.0280S5 

0.028055 
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. 2 
m. 

0.0088247 

0.00860 

. 2 
in. 

0.013953 

0.013953 

lbm/sec 

0.220 

0.220 

lbm/sec 

0.003663 

0.003663 

- 

0 

0 

lbm/sec 

0 
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16.86 

in. 

5.12 

in. 2 

5.0 

in* 


Full Taper 

(Area change linear with distance) 
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SAMFLE CASE ELEMENT AND CI1AMGER GEOMETRY 
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Chamber Geometry 
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